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CHAPTER  1 


INTRODUCTION 


The  nonlinear  interaction  of  waves  in  plasma  has  been 
studied  extensively  Cl-9].  Recently,  it  has  been  shown  that 
certain  nonlinear  wave-wave  interacting  systems  with  a  small 
number  of  degrees  of  freedom  exhibit  chaotic  behaviors  [10- 
153.  Studies  of  this  phenomenon  may  lead  to  new  theories 
for  turbulence  in  fluids  and  plasmas. 

In  the  classical  Landau-Lif shitz  theory  of  turbulence 
[16],  the  physical  variables  describing  a  turbulent  motion 
are  quasi-periodic  functions  of  time: 

x(t)  =  f(u»1t,  u»2t,...,  u>Nt),  (1.1) 

where  f  is  composed  of  periodic  motions  with  frequencies 
u)^,...,  which  are  irrationally  related.  As  the  number  N 
of  modes  becomes  large,  x(t)  generally  exhibits  a  very  com¬ 
plicated  behavior.  The  motion  described  by  (1.1)  can  be 
produced  by  a  sequence  of  successive  bifurcations  of  peri¬ 
odic  motions.  As  certain  parameters  vary,  a  stable  equilib¬ 
rium  point  becomes  unstable  and  a  stable  periodic  orbit 
appears.  Then,  the  periodic  orbit  becomes  unstable  and  a 
stable  periodic  orbit  with  two  frequencies  appears.  This 
bifurcation  process  continues  on. 

It  is  known  that  a  quasi-periodic  function  of  time 
does  not  have  the  mixing  property  (i.e.,  the  autocovariance 
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R ( t )  of  x  does  not  decay  to  zero  as  t  -*•  <*>).  This  implies 
intuitively  that  the  temporal  variation  of  x  is  net  sensi¬ 
tive  to  initial  conditions,  thus  contradicting  an  important 
physical  feature  of  turbulent  states.  Furthermore,  a  quasi- 
periodic  function  is  unstable  in  the  sense  that  a  small  per¬ 
turbation  of  the  parameter  could  lead  to  a  periodic  motion 
[17].  Hence,  the  quasi-periodic  motion  does  not  provide  a 
satisfactory  description  of  the  turbulent  states.  It  is  also 
known  that  some  physical  systems  exhibit  abrupt  transition 
from  a  nonturbulent  state  into  a  turbulent  state  without  the 
appearance  of  periodic  motions.  Such  phenomenon  cannot  be 
explained  by  the  above  mentioned  mechanism. 

Recently,  Ruelle  and  Takens  [18]  proposed  the  following 
definition  for  a  turbulent  state:  In  this  state,  the  trajec¬ 
tories  of  the  system  model  are  attracted  to  a  nonempty  set 
called  the  "strange  attracter"  which  is  a  positive  limit  set 
of  trajectories  consisting  of  neither  periodic  orbits  nor 
equilibrium  point.  On  this  attractor,  the  system  trajec¬ 
tories  behave  in  a  chaotic  manner  resembling  a  turbulent 
state . 

According  to  the  above  notion  of  a  turbulent  state,  it 
is  possible  that  simple  models  described  by  finite  dimen¬ 
sional  systems  of  ordinary  differential  equations  could  have 
chaotic  behavior.  A  well  known  example  is  the  Lorenz  model 
[19]  describing  the  convection  in  a  fluid  layer  forced  by  a 
linear  temperature  gradient: 
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(1.2) 


=  -  ax^  +  ax2 , 
x0  =  rx^  -  x2  -  x2x3  ’ 

x3  =  “  bx3  +  X1X2 ’ 

where  a,  r  and  b  are  constant  real  parameters.  The  Lorenz 
system  exhibits  abrupt  transition  into  a  very  complicated 
behavior  when  the  parameters  pass  through  certain  threshold 
values.  Many  r  :merical  experiments  have  been  performed  on 
this  system  and  their  results  suggest  the  existence  of 
strange  attractors  and  the  mixing  property,  but  they  have 
not  been  verified  mathematically.  In  this  work,  the  term 
"chaotic  solutions"  refers  to  those  solutions  which  lie  on  a 
strange  attractor.  If  a  solution  exhibits  complicated  behav¬ 
ior  but  without  a  mathematical  proof  as  being  a  "chaotic 
solution",  it  will  be  referred  to  as  a  "pseudo-chaotic  solu¬ 
tion".  Hence,  "pseudo-chaotic  solutions"  may  be  periodic 
solutions  with  long  periods,  or  may  correspond  to  trajec¬ 
tories  converging  to  a  stable  closed  orbit  or  an  equilibrium 
point  in  a  complicated  manner.  Many  simple  mathematical 
models  having  chaotic  or  pseudo-chaotic  solutions  have  been 
found  in  such  areas  as  geophysics  [20],  biology  [21,22], 
chemistry  [23]  and  laser  physics  [24]. 

1.1  Turbulence  in  Plasmas 

In  the  conventional  theory  of  turbulence  in  plasmas, 
turbulence  is  described  as  the  state  in  which  a  large  number 
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of  collective  degrees  of  freedom  are  strongly  excited  [25]. 
This  means  that  the  energy  of  the  collective  modes  is  much 
larger  than  the  thermal  fluctuation  level,  and  that  the 
number  of  modes  is  so  large  that  the  plasma  has  a  complicated 
behavior.  The  energy  of  the  unstable  modes  is  distributed 
to  other  modes  by  nonlinear  processes  and  then  dissipated 
through  some  form  of  damping  mechanisms.  When  the  energy 
transfer  between  the  modes  is  balanced,  a  stationary  broad 
spectrum  corresponding  to  stationary  turbulence  appears. 

When  the  energy  of  each  mode  is  small  and  the  nonlinear 
coupling  between  various  modes  is  weak,  and  when  the  linear 
growth  and  damping  terms  can  be  neglected,  the  random-phase 
approximation  can  be  used  to  obtain  the  turbulent  state. 

In  this  approximation,  it  is  assumed  that  the  phase  of  each 
mode  varies  randomly.  On  the  other  hand,  if  the  linear 
growth  or  damping  terms  are  not  negligible,  the  phase  of  each 
mode  is  no  longer  independent  even  when  the  nonlinear  coup¬ 
ling  of  modes  is  weak.  In  spite  of  this  limitation,  the 
random-phase  approximation  is  often  used  for  expedience. 

As  mentioned  earlier,  the  existence  of  a  large  number 
of  modes  does  not  necessarily  mean  that  the  system  is  in  a 
turbulent  state.  Furthermore,  the  assumption  of  independent 
phases  between  the  modes  is  usually  made  without  justifica¬ 
tion.  At  the  initial  stage  of  turbulent  states,  it  is 
reasonable  to  assume  independent  phases  since  the  randomness 
of  thermal  fluctuation  at  the  initial  time  is  retained.  On 


the  other  hand,  once  the  energy  of  turbulence  is  far  from 
the  equilibrium,  there  must  be  some  mechanisms  which  maintain 
the  random  phases.  Flynn  and  Manheimer  [10]  showed  that 
even  when  the  modes  have  no  linear  growth  or  damping  over  an 
interval  of  wave  numbers,  the  phases  are  completely  coupled 
when  the  linearly  growing  and  damping  modes  are  introduced 
outside  the  interval.  This  means  that  the  random-phase 
approximation  is  not  suitable  even  in  the  interval  over  which 
both  the  linear  growth  and  damping  are  negligible. 

At  present,  there  are  a  few  known  models  in  plasma 
physics  having  a  small  number  of  modes  which  exhibit  chaotic 
behaviors  [11-15].  Flynn  and  Manheimer's  model  [10]  consist¬ 
ing  of  ten  modes,  two  of  which  are  linearly  growing  and 
damped,  has  chaotic  behavior.  Rabinovich  and  his  coworkers 
introduced  following  two  simple  models  of  wave-wave  interac¬ 
tion  having  chaotic  motions  [11-13]. 

A1  -  Y-^  =  -  i{A2  +  A1(v1|A1|2  +  w1|A2|2)}, 

•  (1.3) 

A2  +  Y2A2  =  "  ^A1A2  +  A2  ( v2  |  Ai  |  ^  +  w2  I A2  I  ^  ^  » 

and 

A^  +  =  “  i(A2A3  +  vAqA*), 

A2  +  Y2A2  =  “  ^^A1A3  +  v^oAl^’  (1.4) 

A3  +  y3A3  =  -  iA^A* , 

where  A^  is  a  normalized  complex  wave  amplitude;  all 
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parameters  are  real  numbers  and  y^O.  System  (1.3) 
describes  the  interaction  of  two  plasma  waves  belonging  to 
different  branches.  Chaotic  motion  appears  only  when  a 
small  frequency  mismatching  is  introduced.  System  (1.4) 
describes  a  three-wave  interaction  such  that  an  external 
wave  with  a  constant  amplitude  decays  into  two  waves  which 
have  another  resonant  wave.  Numerical  experiments  for  these 
systems  suggest  that  they  have  strange  attractors  and  mixing 
properties.  Hence,  they  could  represent  turbulent  states. 

We  should  note  that  in  these  systems,  the  linear  instability 
may  lead  directly  to  a  turbulent  state,  and  that  the  chaotic 
behavior  is  intrinsic  to  the  equations  and  no  additional 
assumption  of  randomness  is  necessary. 

1.2  Outline  of  Dissertation 

In  this  dissertation,  three-wave  interacting  systems 
having  chaotic  behaviors  are  studied  in  detail.  Asymptotic 
phase  locking  of  the  waves  as  t  -*■ 00  are  also  studied,  which 
is  a  special  property  of  wave-wave  interaction  with  linear 
damping  terms.  Once  the  phase  are  locked,  the  system 
reduces  to  a  three-dimensional  one. 

In  Chapter  2,  the  systems  to  be  studied  in  this  work 
are  presented.  All  these  systems  are  obtained  from  the 
equations  for  a  set  of  interacting  oscillators.  The  reduced 
three-dimensional  equation  is  obtained  for  each  system.  We 
also  give  sufficient  conditions  for  the  reduced  equations  to 
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describe  the  asymptotic  behavior  of  the  original  system. 


In  Chapter  3,  the  reduced  equations  are  studied  mainly 
by  a  root  locus  analysis  of  the  characteristic  equations  of 
the  linearized  vector  field  at  the  equilibrium  points.  From 
the  results,  we  choose  those  systems  whose  reduced  equations 
may  have  periodic  or  pseudo-chaotic  solutions  for  almost  all 
initial  conditions.  Based  on  physical  considerations,  two 
systems  are  selected  for  detailed  study. 

In  Chapter  4,  the  results  of  numerical  experiments  for 
two  chosen  systems  are  presented.  The  transition  between 
periodic  and  pseudo-chaotic  motions  are  shown  in  detail. 
Simple  two-dimensional  models  are  introduced  which  are  used 
to  explain  the  transition  of  behaviors. 

In  Chapter  5,  we  obtain  the  first  return  mappings  of 
the  trajectories  for  the  chosen  systems,  some  of  which 
resemble  one-dimensional  mappings.  We  discuss  certain 
statistical  properties  of  general  one-dimensional  mappings 
and  the  first  return  mappings  of  our  systems. 
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CHAPTER  2 


WAVE-WAVE  INTERACTIONS 


In  this  chapter,  we  introduce  the  systems  of  wave-wave 
interactions  to  be  studied  in  subsequent  chapters.  All 
these  systems  are  derived  from  the  equations  describing  a 
set  of  nonlinearly  interacting  oscillators.  A  wave-wave 
interaction  which  does  not  belong  to  this  framework  is  not 
considered  here. 

2.1  Wave-Wave  Interacting  Systems  Y =  F(Y) 


We  begin  with  the  equations  for  a  set  of  nonlinearly 
interacting  oscillators  described  by 


V.,  x.x. 


+  I  viklxixkxl’  ^ 6  ^N ^  {!>••• »N)» 

j,k,leZN  3kl  3  k  1  N 


(2.1) 


where  x^  is  the  displacement,  uk  is  the  natural  frequency, 
and  is  the  growth  or  damping  coefficient  for  the  i-th 
oscillator  ;  v.,  and  v.,  ..  are  the  nonlinear  coupling 

]  K  ] 

coefficients  of  the  oscillators.  All  the  variables  and 

parameters  are  real  numbers.  We  assume  that  the  linear 

damping  and  growth  rates  are  small  (i.e.,  <<  io^)  and  the 

nonlinear  terms  are  also  small  (i.e.,  |v.,x.|,  |v.,,x.x,  I 

]k  ]  jkl  3  k 1 

2 

<<  ok  ),  Then,  we  can  express  x^(t)  in  the  form: 


(2.2) 


x.(t)  =  a.  (t)exp(ico.t)  +  a*  (t)exp(-iw. t)  , 

1  1  l  l  l 

1£ZN’ 

where  a^  is  a  slowly  time-varying  component.  We  also  assume 
frequency  matching: 


co  .  =  E  U  .  co . 
1  icZ  3  3 


Pj  e  (-1,1} ,  i «  ZN, 


(2.3) 


where  uk  >  0  and  uk  t  ok  if  i  i  j  .  Then,  the  frequency  matching 
conditions  for  N=  2,3,4  are  shown  in  Table  2.1,  where  it  is 
assumed  that  the  v...  's  in  (2.1)  are  zero  for  N=  3,4.  For 
N  =  3,  the  first  frequency  matching  condition  has  two  degrees 
of  freedom  and  the  remaining  ones  have  only  one  degree  of 
freedom.  For  N  =  4,  we  assume  that  one  of  the  waves  is  an 
external  wave  with  fixed  amplitude.  Then,  there  is  one 
degree  of  freedom  in  the  first  three  cases  and  none-  in  the 
rest.  The  energy  of  the  waves  is  obtained  by  an  integration 
over  the  frequency  domains.  Here,  for  N=  3,  the  energy 
contribution  of  the  last  two  cases  is  negligible  as  compared 
to  that  of  the  first  case.  In  what  follows,  we  consider 
only  the  first  case.  Similarly,  for  N=  4,  only  the  first 
three  cases  will  be  studied. 

For  N=2,  substituting  (2.2)  into  (2.1)  and  normalizing 
the  variables,  we  have 


A.  +  y-A.  =  -i{vA .  2  +  A .  ( V  .  .  .  I A .  |2  +  V...IA.I2)},  } 
i'll  3  1  111 1  1 1  13  3  3  ! 

A.  +  y.A.  =  -i{A.A*  +  A. (V. .  .  |A.  | 2  +  V.  .  .  I  A.  I  2)}  , 

3  3  3  13  3  3H1  11  333  1  3  1 


(2.4) 
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Table  2.1  Frequency  Matching 


KB 

(o).,w.)  =  (  2  a) .  ,  ii> .  ) 

1  3  3  3 

N  =  3 

,uik) 

(u,j+tl,k’“j  ,uk> ,  (2ok  ,3o>^  )  ,  (4wk,2aik,4)k) 

N  =  4 

(  ou)k+vai^  ,uj^/  2  ,uik  )  , 

(2oj^  ,C4Jk+V(u1  ,a)k,4J1> ,  (o,v)  e  {  ( 1 , 1 )  ,  ( 1 , -1 )  ,  ( -1 , 1 )  }  , 

(u)j+ajkla)k+oai1,a»k,u1)  ,  ae{-l,l}, 

(8u)1,4<o1>2cj1,<»)1)  ,  ( 4aik ,  2tok ,  wk ,  3wk ) , 

(4a)k,2tok,a)k,6u>k)  ,  ( 4u)k  ,  2uh.  ,u>k ,  5^)  , 

(2oj.  ,aj.,3oj.,(3/2)aj.).  (2cu.,cj.,6aj.,3{j.) 

3  3  3  3  3  3  3  3 
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where 


ic 

A  .  =  ( v .  .  /  2uj  .  )a  .  ,  A.  =v(  I  v  .  .  v  .  .  I /2w  .a) .  )^a  .  ,  v  =  sgn(  v  .  .  v  .  .  )  , 
i  13  113  1  33  13 1  13  3  6  33  13 

V.  .  .  =  {  (  2u> .  /  v  .  .  )  2/2o).  }v.  .  .  ,  V.  .  .  =  {(4u>.w./|v.  .  v .  .  |  )  /  2  w .  }  v .  .  .  , 
111  3  3  3  1  in  133  13'  13  3  3  1  133 

V.  .  .  =  {  (  2u> .  /  v .  .  )  2  /  2ui .  }v.  .  .  ,  V.  .  .  =  {(4a).d)./jv.  .v.  .  |  )  /2u> .  }  v  .  .  .  . 
311  3  13  3  3ii  333  1  3  1  13  33 1  3  333 


For  N  =  3  , 

A. 

we  have 

+  Y  .  A .  = 

-  iv.A.A. 

1 

'  1  1 

1  3  k 

A. 

+  Y -A .  = 

-  iv.A*A. 

3 

3  3 

3  k  1 

•  < 

+  y,  A.  = 
'k  k 

-  iv.A.A* 
k  1  : 

where 

Ai=(ivij',kii/4“j“k)’!ai’  Ar(ivijvjki/4“i“k)'V 

V<lvkivjkl/“"i“j)'V  vj=sen<vijuki)- 

For  N  =  4 ,  we  have  three  cases  as  follows: 

Case  (i) : 


A. 

+ 

Y  -  A  .  = 

-  i  ( v !  vA .  A.  ) , 

% 

1 

1  1 

1  3k’ 

> 

A. 

+ 

y.a.  = 

-  i ( v ! vA . A*  + 

v  .A,  A.  ) 

3 

T  3  3 

3  ik 

3  k  1 

<r 

>  • 

+ 

YkAk  = 

-  i(v'vA.A$  + 
K  13 

V:A*> 

Ai  *  YiAi 


- 


where 
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(2.5) 


(2.6) 


J 


v<|vijvik|/,“j“k),Vv>  Aj=<|vjkvjii/4“k“i)S’ 

Ak=<lvkivjkl/'*“j"i>\’  Ai  =  (lv5kvjil/'*"j“k)S’ 

v=<lvji/vijl“i/“i>l5’  vrS|?n(vijvik)*  ''rssn(vijvjk)’ 

uk=sgn^vikvjk> ’  v j =sgn< vjkvH 1 ’  uk=s«"(vklvjk)- 

vl=sgn(vjlvkl). 

Here,  without  loss  of  generality,  we  set  w .  =  w .  + w,  ,  w .  = 

1  3  k  3 

w,  +  to,  ,  since  if  o  =  -l,  the  substitutions  A.-*- -A.,  A.-*- -A., 

k  1  i  i  ’  j  ]  ’ 

A^  ■+■  -A*  and  -*■  lead  to  the  same  equation  as  (2.6). 

Case  (ii)  : 


A.  +  y . A. 
i  'i  l 


=  -  i ( v ! vA?  +  v -A  A, ) , 
11  i  k  1 


A .  +  y .A. 
1  3  3 


\  4  YkAk 


iCvAjA*) , 


=  -  Kv^A*), 


(2.7) 


A1  +  Y1A1 


i(vlAiAk)» 


where 


Ai=<|vilvikl/““k“l>  ai’  Aj=(|v..v..|/«„i»j)'5aj/v, 

V<  1  vikvkl  I  ^»1»1>'5aK  -  *i  =  <  I  vilvkil - 


v=<vij/2“j)<4“kVlVi1vikJ^.  Vsgn(vikvil)’ 


v;  =  sgn<v  v  >,  vk=sgn(v.kvkl), 


and  a)  .  =  2ui .  =  cj,  +  u  . 
i  ]  K  l 
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Case  ( iii) : 


A.  +  y .A. 
1  11 


A.  +  y.A. 
3  3  3 


A,  +  y.  A, 
k  k  k 


A1  +  Y1A1 


where 


i( v . vA? )  , 
i  3 

iCvA.A*  +  v  .A.  A,  )  , 
13  3  k  1 

1<VjA‘)> 

i(ViA.A*>, 


(2.8) 


A.  =  (vij/2wj)aj/v,  A.  =  (|vjkvjl|/4a)kuJl)55aj, 

Ak=(lVjkVkll/4(iljwl)l5ak’  Al=(lvjlvkll/4u)jwk)l5al’ 
v'(  lvijv5j  l/4“i“j>  <,*'V,1/Ivjkvnl)  >  vr-sgn(v..vjj) 
vj  =  sgn(Vj^Vji) ,  vk=sgn(v.kvkl),  v^sgnCv.^  > , 


and  u> .  =  2u) .  ,  w .  =  u,  +  u),  . 
l  3  3  k  1 

As  mentioned  earlier,  we  fix  the  amplitude  of  one  of 
the  waves  as  an  external  wave  for  N=  4.  We  choose  the 
external  wave  so  that  the  equation  is  nonlinear,  since  linear 
equations  do  not  have  chaotic  solutions.  Then,  from  (2.6), 
(2.7)  and  (2.8),  we  obtain  five  distinct  three-wave  interac¬ 
ting  systems.  Thus,  we  have  seven  systems  which  may  have 
chaotic  solutions  and  they  are  listed  in  Table  2.2.  Here, 
without  loss  of  generality,  A^  is  taken  to  be  real.  Also, 
the  substitutions  An-*-An/v,  A  -*■  A  /v,  m=i,j,l  and  v-*-l/v 
have  been  used  to  derive  (51)  and  (7'). 

In  Table  2.2,  some  systems  with  certain  values  of 
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Table  2.2  Wave-wave  Interacting  Systems  Y  =  F(Y) 


(wi  ,UK  ) 

=  (  2uj  .  ,  a) .  ) 

3  3 

A.+Y-A.  = -i{vA.2+A. (V. . . |A. I 2+V . . . |A. I2)}, 
ill  3  l  in 1  i 1  ID  3  ] 

A.+y.A.  = -i{A.A*+A.(V. . . |A. |2+V. . . |A. |2)}, 

32D  133  311  11  333 1  31  ’ 

V  £  {-1,1}  (1») 

(60i,0)j  ,0)k> 

=  (u,j+tV 

wj’uk) 

Ai+YiAi =  "iviA3Ak’  vm  6  {-1,1} 

A.+y.A.  =  -iv.A.A*, 

333  31k’ 

VW-iViA!  <2'> 

(  u>i  ,  Ui .  ,0)k, 

0)1)  =  (dK+O^, 

(jJk+a)l,aJk’0Jl) 

A:*yjAj  =  -1(vAVvjvVk)-  Ao  *  Ai« 
VVk  2  -i<vkAjAi*',k,,AoAj)>  V  f-1-11’ 

AlnlAl  =  -ilviAjAk>  <3,) 

o>i )  =  («k+«1» 

uk/2  ,uk,0J1) 

VYiAi  ’  -i(VkV',i'V>’  Ao  =  Aj’ 
VYkAk  =  -i<ViAl>> 

A1*Y1A1  '  '“’M1 

VYiAi  =  A0  =  Ak, 

Aj+YjAj  =  -i<AiA?),  .,<{-1,1), 

A1  +  y1A1  =  -Hv^vAqA^  (5’) 

(aJi’UJj’Wk 

ai-^ )  =  (  2o)  j  , 

<Jk+a>l’a)k’a,l) 

yY5Ajs-i<vAV''v?)>  Ao=Ai’ 
AktYkAk  =  -i(vkAjAi>>  v1'1’11’ 

A1+Y1A1 = ■i(vlAjAk)  (6’> 

A.+y.A.  = -i(v!A.2) ,  A  =  A,  , 

111  1  ^  U  iC 

A^+YjA.  =  -iCA^+Vj vA[)A1)  ,  vme  {-1,1}, 

ai+yiai  =  -i(vivAoAj)  (7>) 
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parameters  obviously  have  no  chaotic  solutions.  For  example 
consider  system  (3')  with  y .  , y.  >  0 ,  y,  <  0  ,  v  =1,  m  =  j ,k , 1 
and  v^=  1,  m=  j,k.  By  direct  computation,  we  have 

(d/dt ) ( |A. |2-|Ak|2+2|A1|2)/2 

=  -Yj |A. |2-Yk|Ak|2-2y1|A1|2.  (2.9) 

Hence,  the  quantity  ) A^ ( t ) | 2 - j Ak ( t ) | 2  +  2 ) A^ ( t ) | 2  decreases  as 
t  increases  and  therefore,  the  trajectory  of  (3')  diverges 
to  infinity  as  t  -*■ 00  for  almost  all  initial  conditions.  In  a 
similar  manner,  we  can  discard  certain  systems  in  Table  2.2, 
and  only  those  cases  in  Table  2 . 3  may  possess  chaotic 
solutions . 

In  what  follows,  we  use  subscript  1,2,3,...  in  place 
of  i,j,k,...  .  For  example,  subscripts  j ,k  and  1  of  (3') 
are  replaced  by  1,2  and  3,  respectively.  Moreover,  we 
denote  systems  (2')  -  (7')  by 

Y  =  F(Y),  (2.10) 

A  T 

where  Y  =  ^yi>zi>y2’z2’^3’z3^’  and  Yj  =  Re (A^ )  and  z  ^  =  Im(A^  ) 
j  =  1,2,3.  We  also  denote  the  solution  of  (2.10)  at  time  t 
corresponding  to  initial  condition  YQ  at  t  =  0  by  '/(t,YQ), 
and  the  positive  semi-orbit  by  /+(Yq)  =  {Y  :  Y  =  Y(t,YQ)  ,  t  2  0} 

2.2  Reduced  Systems  X =  F(X) 

System  (!')  in  Table  2.3  is  identical  to  (1.3).  That 
is,  model  (1.3)  [11]  is  the  only  two-wave  interacting  system 
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Table  2.3  Parameters  for  Chaotic  Solutions 


( 1 '  ) 

v  =  1 

(sgn(yi) ,sgn(Yj ))  =  (-1,1) 

(2  '  ) 

(vi’vj’V 

(sgn(Yi) , sgn( Yj ) ,sgn(YR) ) 

(1,1,1) 

(-1,1,1) 

(3'  ) 

( vj  ,  v£  )  ,  (  V  j  ,  vk ,  v1 ) 

( sgn ( Yj ) , s  gn ( Yk ) ,sgn(Yl)) 

(1,1),  (i,i,i) 

(±1,±1,±1),(1,±1,-1),(-1,±1,1) 

(1,-1,-1) 

(1,1, ±1) 

(-1,1,-1) 

(1,1, ±1) 

(-1,-1,1) 

(±1,±1,±1) ,(±1, 1,-1) , (±1,-1, 1) 

(1,-1),  (1,1,1) 

(1,-1, ±1) ,(-1,1, ±1) 

(1,-1,-1) 

(±1,±1,±1) ,(1,-1, ±1) ,(-1,1, ±1) 

C-1,1,-1) 

(1, ±1,-1) ,(-1,1, ±1)  ,(±1,-1,1) 

(-1,-1, 1) 

(1,-1, ±1) ,(-1,1, ±1) 

(5'  ) 

(v[,vi,v1) 

(sgn(Yi) , sgn ( y j ) ,sgn(Y1) ) 

(1,1,1) 

(1,±1,-1),(-1,1,±1),(±1,-1,1) 

(1,1,-1) 

(±1,±1,±1) ,(1,-1, ±1) , (-1,1, ±1) 

(-1,1,1) 

(-1,-1,-1) ,(1, ±1,-1) ,(-l, ±1,1) 

(-1,1,-1) 

(±1,±1,±1) ,(±1,1,-1) ,(±1,-1,1) 

(7*  ) 

(  v  !  ,  v  .  ,  v ,  ) 

( sgn ( y^ ) ,sgn(Yj ) ,sgn(Y1) ) 

(1,1,1) 

(1,-1,±1) ,(±1,1,-1) ,(-l,±l,l) 

(1,1,-1) 

(±1,±1,±1),(1,-1,±1),(-1,1,±1) 

(-1,1,1) 

(-1,-1,-1) ,(1,±1,-1) ,(-l,±l,l) 

(-1,1,-1) 

(±1,±1,±1) ,(±1, 1,-1) ,(±1,-1, 1) 
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which  may  have  chaotic  solutions.  Equation  (1')  is 
rewritten  as  a  three-dimensional  equation: 


^  =  -y1x1-6(ux2+2x1x2-x2{(V111-2V211)(x12+x22)  j 


+(V122~2V222)x3*’ 

x2  =  -y1x2+<Swx1+x3-2x12+x1{(V111-2V211)(x12+x22) 


+(V122~2V222^X3^’ 


(2.11) 


X3  =  -2y2x3~2x2x3,  x3>0,  J 

2 

where  =  r^os^f ,  x2  =  r^inffr,  x3  =  r2  ,  r1  =  jA-J  ,  r2  =  |A2  |  , 

(J>  =  2Arg(A2  )-Arg(A^ )  and  6to  is  a  frequency  mismatching  u)^-2a)2. 
Equation  (2.11)  has  pseudo-chaotic  solutions  with  or  without 
the  third  order  terms,  and  its  behavior  has  been  studied 


[11,13,14]. 


Now,  we  shall  show  that  systems  (2')- (7')  also  reduce 
to  three-dimensional  equations  on  certain  subsets  of  the 
original  six-dimensional  space.  We  choose  equation  (3')  as 
an  example .  Let  A j  =  r . exp ( i0  ^ ) ,  r  ^  >  0  ,  j  =  1,2,3.  Then , 
equation  (3')  can  be  written  as 


r1  +  ir101+Y1r1  =  -itv^r^xp  {-i ( 0 1"02-0  3  ) 

+vivAor2exp{-i(61  +  02)}  J, 
r2+ir202+Y2r2  =  -i[v2r1r3exp{i(01-02~03 ) 
+v2vA3r^exp(-i(0^+02)}], 


(2.12) 
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Vir3*3+Y3r3  =  -i^3r1r2exp{i(  8  1-02-0  3  )}.  J 


Let  us  define  the  sets: 


£q  =  {Y  :  cos ( 0 ,-0 2-0 3 )  =  cos(03+@2)  =0,  0^  =  0, 


rj  *  0,  j  =  1,2,3}, 


£  =  uL 

r  0  9 


£n={Y:r^  =  0,  i  =  1 , 2  and/  or  3  }  , 


£0  =  :  ri  =  r2  =  0}, 


£0  =  {Y:r^  =  r.=0,  i =  1,2  or  3}  -  £Q, 


Z0  =  Zn  "  *Z0  u  Z0*  ‘ 


Obviously,  £Q  is  not  closed,  and  £  and  £  u £Q  are  closed. 

0  n  n  0 

Furthermore,  £  n  £  =  <p .  Hence,  £Q  -  £Q  is  nonempty  and 

n  0  u  w 

£0  -  £q  c £  An  example  of  £0  is  shown  in  Figure  2.1.  The 

set  £^  is  located  at  the  origins. 

We  assume  that  Y(0,Yn)=Yne£.  .  From  (2.12), 

u  u  y  o 

0.(t)  =  0,  j  =  1,2,3,  if  y(t,Y  )  €  £  .  Hence,  Y+(tn)  c  £„  or 

J  U  U  o  U  t)  o 

there  exists  a  t.  >0  such  that  Y(t,Yn)e£Q  for  t  «  [ 0 , t, ) 

1  U  U  o  1 

and  y(t,  ,Yn)  €  £.  -  £  <=  Z  .  Since  £  =£no£  u£  and  £n  is 

xu  u  o  ”  o  **  nuuu  u 

an  invariant  set  (i.e.,  if  YQ  e  £Q ,  Y( t , Yq )  £  £Q  for  t  >  0), 
y(t3,Y0)  £  £0  uE0.  Suppose  that  /(t1,Y0)  £  £0.  Then,  YCt-^jY) 
is  located  at  one  of  the  three  origins  in  Figure  2.1.  Since 
F(Y)  is  transverse  to  £0  at  Y =  Y(t^,Yg)  and  F(Y)  of  (3')  is 
differentiable,  Y (t.+e,Y_)  e  Za  for  arbitrarily  small  e  >  0. 

1  U  U  0 
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Figure  2.1  Phase-Locked  Trajectory 


Suppose  that  y(t^,YQ)e  £q.  Then,  y(t^,Yg)  is  located  at  two 
of  the  three  origins  in  Figure  2.1.  For  F(Y)  of  (3'),  there 
is  no  t>0  such  that  y(t,YQ)e£Q  for  te  [t-^t  +x).  Moreover, 
F(Y)  of  (3')  is  continuously  differentiable.  Hence, 
y(t,+e,Yn)  e  £.  for  arbitrary  small  e  >  0.  Thus,  the  set  £„ 

i  u  y  o  y  o 

is  an  invariant  set,  and  the  set  Zn  is  invariant  in  the 

y  o 

sense  that  if  Yn  e  I  ,  y(t,Yn)  e  Za  for  almost  all  t  except 
u  y  o  u  y  o 

when  r.(t)=0,  i  =  1,2  and/or  3 .  Thus,  V+(Yn),  Yn  e  £Q  ,  is 
l  u  u  o  o 

contained  in  £«  - 
y  o 

•»  A  — 

If  y(0,YQ)  e  (Zq  u£q)  n  a^ove  used  arguments, 

y(e,Yn)  «  £.  for  arbitrarily  small  e  >  0.  Here, 
u  y  o 


9n  =  lim  9  ,(t)  . 
U  t-*-0 


Such  a  limit  0q  exists,  since 
Hence,  y+(Yn)  c £  if  Yn  e  (£n 


a  set  £ 


r0  ' 


£  _  =  £  u {(£n  u  £n)  n  £  } 
r0  r  0  0  r 


F(Y)  of  (3')  is 

u  £n )  n  £  .  Thus 
0  r 


5 


continuous . 

,  if  we  define 
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y  +  ( Yn  )  c  E  for  all  Yn  e  I!  n ,  where  if  Yne  I  ,  0  n  =  0 .  (  0  )  and 
0  0o  0  ru  ’  Or’Ol 

if  Y0  £  d0  U  V  °  V  9 0  =  iimn  e3(t)- 

t-*-  0 

Now,  we  introduce  a  reduced  equation  of  (3').  If 

y(°’Yo>^0o> 

0i(t)  =  0^0)  +  imi(t),  i=  1,2,3,  (2.13) 

where  63(0)  =  6g  and  n^  is  an  integer-valued  step  function 

having  discontinuities  only  when  r^(t)  =  0.  Let  us  define  a 

A  T 

vector  X  Q  (x.»x2,x3)  : 


x^(t)  =  var1(t)exp{in1(t) } , 
x2(t)  =  vr2(t)exp{in2(t) } , 
x^(t)  =  var3 (t)exp{in3 (t) } , 


(2.14) 


where  v  =  sin( 0.^(  0 )-02  ( O)-0 3 ( 0 ) )  and  a  =  sin(  0.^ ( 0 )  +©2  ( 0 ) ) . 
Then,  x.  is  real,  and  at  the  time  t  when  Y(t,Yn)  e  £  , 

1  u  b  0 

equation  (2.12)  reduces  to 


x2  =  -  Y1X1  -  v1x2X3  "  vivA0X2’ 
x2  =  -  Y2x2  +  V2X1X2  "  v2vA0X1’ 
x3  =  “  Y3X3  +  V3X1X2- 


(2.15) 


This  equation  actually  holds  at  the  time  t  when  y(t,Yg)  « 

Eq  -  la  ,  since  its  trajectory  stays  in  Ea  -  £Q  momentarily. 
00  w  0  0000 

If  Y.  £  !  n  (En  uZn),  the  solution  of  (2.15)  is 
Or  0  0 

identical  to  that  of  (2.12).  For,  the  trajectory  of  (2.15) 
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momentarily  stays  in  the  set  in  the  X-space  corresponding  to 

_  /v  ~ 

n  (Zq  uEq).  Furthermore,  equation  (2.15)  has  a  solution 

X(t)  = 0V,  the  zero-vector  in  the  X-space,  or  X(t)  =  (0,0, 

X 

T 

x3 ( 0 )exp(-'f  gt )  )  on  the  set  in  the  X-space  corresponding  to 
E  ,  and  they  are  identical  to  those  of  (2.12)  on  Eg.  Thus, 
equation  (2.12)  reduces  to  (2.15)  for  all  Yg e  E^g  uEg.  ^ 

Y0  e  lr0  -  V  ^  V  ■=  Sr0  *  Eo  and  if  Yo  ‘  V  <=  V 

Similarly,  each  of  the  remaining  systems  in  Table  2.2 
reduces  to  a  three-dimensional  equation  on  a  certain 
invariant  set.  These  equations  are  given  in  Table  2.4  along 
with  the  definitions  of  variables  and  E  .  Systems  (2)-(7) 
and  their  solution  will  be  denoted  by 


X  =  F(X),  (2.16) 

and  X(t,XQ),  respectively,  where  XQ  is  an  initial  value. 
Moreover,  we  denote  the  positive  semi-orbit  by  X+(Xg) 

=  (X  :  X  =  X(t,XQ) ,  t  >  0}. 

For  some  reduced  equations,  fx(t)|[  <*»  or  0  as  t  •+ 00  for 

initial  points  in  certain  open  subsets  of  the  X-space.  Let 
us  consider  (3)  with  (v^,v3)  =  (1,1),  (vi>v2’v3^  =  (1,1,1)  and 
( sgn( Yji  ) » sgn(  y2  )  ,  sgn( y3  ) )  =  (1,1, -1)  as  an  example  (see  Table 
2.3).  It  is  not  determined  whether  the  original  system  (3') 
with  such  parameters  has  a  diverging  trajectory.  On  the 
other  hand,  from  (3), 

(d/dt){(x3-(3/2)vAg)^-x1^-2x22} 
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sin(ei(0)-2eo(0)) 


=  Y,  x 
1  x 


'  +2y  2x2 


2+Y3(x3-(3/2)vA0)2-(9/4)y3 


vA, 


2  2 

Hence,  if  (0,Xg)  is  inside  the  elliptic  cone  C  :  x^  +2x^ 

2  2  2 
=  (x3- (  3/ 2  )  vAg  )  and  outside  the  ellipsoid  E  :  Y-j^  +2y2X2 

+Y3  (x3-(  3/2  )  vAg  )2  =  ( 9  /  M- )  y  3  I  vAg  |  2  ,  |  X(  t )  ||  -*■ 00  as  t  -*■  00  where 

X(t)  =  (t,Xg)  (Fig.  2.2).  Similarly,  for  some  cases  of  (3) 

and  (4)  in  Table  2.3,  there  exist  open  sets  Rg  such  that  if 

XQeRg,  |  X(  t )  |  ■+■ 00  as  t-+“  where  X(t)  =  X(t  ,Xg  ) .  These 

systems  will  be  discarded  in  Chapter  3,  since  they  do  not 

have  periodic  or  chaotic  solutions  for  almost  all  initial 

conditions.  The  remaining  cases  are  shown  in  Table  2.5. 


x3 


Figure  2.2  Diverging  Trajectory  of  (3) 

2 . 3  Phase  Locking 

In  the  previous  section,  we  have  shown  that  if  the 


phases  are  locked  at  the  initial  time  (i.e.,  y(0,Yn)e  Z  ), 


Table  2.5  Parameters  for  Chaotic  Solutions 


of  Systems  (3)  and  (4) 


■ 

(sgn(y1) ,sgn(y2) 5sgn(Y3) ) 

(i,i)  (i,i,i) 

(1,1,1), (-1,1,1) 

(i,-i, -i) 

(1,1,1) 

(-i,i,-i) 

(1,1,1) ,(1, 1,-1) 

(-i,-i,d 

(1,-1,1) 

(i,-i)  (i,i,i) 

(-1,1,1) 

(-i,i,-i) 

(1,1,-1). 

(-1,-1, i) 

(1,-1,1) 

(4) 

(Vi,V2,V3) 

(sgn(Y1) ,sgn(Y2 ) ,sgn(Y3 ) ) 

(1,1,1) 

(1,1,1) ,(1, 1,-1) 

(-1,-1, 1) 

(1,-1, 1) 

(-1,1,-1) 

(1,1,-1) 

(1,-1,-1) 

(1,1,1) 

¥ 
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they  remain  locked  for  almost  all  t  (i.e.,  Y(t,Yg)e  Eg  for 
almost  all  t  >  0  except  when  Y(t,Y-.)  e.  E  -  -  E  n)  .  Furthermore, 
on  Z^g  uZQ,  the  original  equation  reduces  to  a  three-dimen¬ 
sional  one.  It  is  meaningful  to  study  the  reduced  equation, 
if  it  describes  the  asymptotic  behavior  of  the  original 
equation  in  the  limit  t -*■<*>.  Otherwise,  the  solutions  of  the 
reduced  equation  may  not  resemble  those  of  the  original  one. 
In  order  to  find  whether  the  reduced  equation  represents  the 
original  one  in  the  above  sense,  we  shall  obtain  sufficient 
conditions  for  a  subset  A  of  Z^g  u  ZQ  to  be  an  attractor  in 
the  Y-space.  Here,  we  call  A  an  attractor  in  the  Y-space, 
if  there  is  a  neighborhood  0  of  A  in  the  Y-space  such  that 
if  Yq e  0,  then  y+(YQ)  c  0  and  Y(t)  =  Y(t,YQ)  converges  to  A  as 
t  ”,  that  is , 

lim  inf  ||Y(t)-YS  =  0.  (2.17) 

t-*-°°  YeA 

If  0  is  the  whole  Y-space,  the  phases  become  locked  asympto¬ 
tically  as  t  -*• 00  for  any  initial  conditions.  Even  when  0  is 
a  small  subset  of  the  Y-space,  it  is  theoretically  easy  to 
choose  the  initial  conditions  for  phase  locking,  since  the 
volume  of  0  is  nonzero.  Hence,  if  a  neighborhood  0  exists, 
it  is  meanningfull  to  study  the  behavior  of  the  reduced 
equation. 

We  again  choose  system  (3')  as  an  example.  According 
to  Pikovskii  et  al.  [12],  the  set  E^  is  always  an  attractor 
if  all  of  the  y^’s,  v^'s  and  v|'s  are  positive.  It  can  be 


readily  shown  that  their  analysis  is  inaccurate  (see  Appen¬ 
dix).  Here,  we  shall  give  a  more  detailed  analysis  for  the 
general  cases  of  (3').  Similar  arguments  hold  for  other 
systems  also. 


2.3.1  Differential  Equations  of  Phase  Differences  | =  $£ 

Equation  (2.12)  suggests  studying  the  time  variation 

of  the  quantities  A^A*A*  and  A^A^ ,  since  they  include  the 

phase  differences  anc^  ®l  +  ®2  are  important 

2 

here.  Taking  the  time  derivative  of  A^A*A*,  A^A2  ,  A2A3  and 
2 

A^A* ,  we  have 


’4  J 


■<  VW  0 


-HvAo 


-V2vA0 


-<W 


-2vJvAq  -(2v,|A,|2+v,|A,|2)  -(2y,+Y,)  0 
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'2  3 


[-2v'vAq  -(2v1|A3|2-v3|A1|2)  0 


-(2W. 


or 


|(t)  =  $  ( | A  2  C  t ) |2, ( A  3  C  t ) |  2  )  £  ( t ) , 


(2.18) 

(2.19) 


where  E  ^  =  Re  ( A^*A*  )  ,  ^sReCA^A^,  £3=Im(A2A3)  and 
2 

=  Im(A^A3).  The  system  of  equations  (2.18)  is  not  closed, 
but  it  can  be  regarded  as  a  closed  time-varying  linear 
differential  equation: 


C(t)  =  <Kt)E(t) .  (2.20) 


Let  us  consider  the  set 
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(2.21) 


Z  =  {Y  :  5  =  0  } . 

where  0^  is  the  zero  vector  in  the  £-space.  From  (2.18),  Z 
is  an  invariant  set.  As  shown  later,  if  £(t)-*-0£  as  t -*■“>, 
the  corresponding  Y(t)  converges  to  Z  as  t -*■“>.  Since 
Z  -  (Z  u  Z q )  is  not  an  invariant  set,  Y(t)  does  not  converge 
to  Z-(Z^uZq)  as  t-*-“>.  Hence,  Y(t)  converges  to  Z^  and/or 
to  Zq  as  t  00 .  As  mentioned  in  Section  2.2,  system  (3') 
reduces  to  (3)  on  Z^g  uZ^  and  therefore,  on  Z^  uZ^.  Hence, 
if  £(t)-*-0£  as  t -*■“>,  then  system  (3)  describes  the  asymptotic 
behavior  of  (3'). 

Similarly,  an  equation  in  the  form  (2.19)  is  obtained 
for  each  system  of  (2'),  (4 '  )  -  (  7  '  )  .  Moreover,  if  £  ( t )  -*•  0^ 
as  t  -*■<*>,  each  of  reduced  systems  (2),  (4)  -  (7)  describes  the 
asymptotic  behavior  of  its  original  system.  The  definitions 
for  vector  5  and  matrix  $  for  each  system  are  given  in  Table 
2.6.  Here,  $  of  (2)  is  a  constant  real  number,  and  for  other 
systems,  $  is  a  time-varying  matrix. 

2.3.2  Sufficient  Conditions  for  Convergence  of  £(t)  to  0^ 

For  equation  (2'),  we  see  that  £(t)-+0  as  t  «  if  and 
only  if  y1+Y2+Y3>0>  and  |  £(t)  j  ->■  00  as  t  -►  00  if  and  only  if 
Y^+Y2+Y3<0-  For  the  remaining  cases,  the  following  result 
in  [26]  is  applicable: 

(S.2.1)  Consider  a  time-varying  linear  system 
expressed  by  (2.20)  and  define 
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Is  !!„  =  max  i?il,  U^E®]  = max{$..  + E  |>, 

i  i  j  i  -* 

where  $  =  ($..).  Then , 
i] 

IIS(0)||oc,exp{/J-yco[-$(x)]dT}  <  ||§(t)||co 

<  I!  §  ^  0  >  II  «.exp{ /qU^H  $(x)]dT},  t>0. 

The  above  statement  implies  that  if  u^C^Ct)]  <  0  for 
all  t>0,  then  |  E,  ( t )  ||  w  -*■  0  as  t  -*• 00 .  We  derive  a  modified 
system 

|'(t)  =  $,(t)|’(t)  (2.22) 

from  (2.20),  where  V  =  (p'^,  p^2  ,  P3C3 ,  P^4  )T,  p!  >  0, 
j  =  1,2, 3, 4.  By  applying  (S.2.1)  to  (2.22)  and  since 
»  S'  (t)Ioo  -*■  0  implies  I  S  ( t )  ||  ^  0  ,  H§(t)||oo-^0  as  t if 

yoo[S,(t)]<0  for  all  t>0.  We  consider  system  (3*)  as  an 
example.  For  this  case,  ub[$']  <  0  if  and  only  if 

2v2  lA3 1 2+v3  lA2  I  2  <  (Yl+T2){p3(2Y2+Y3)"2  I  vA0 1  ’ 

1 2vl lA3 1 2_V3 IA1 1 2 1  <  (Yi+Y2){p4(2Yl+Y3)'2lvAo|}/(P3+P4),  ‘  (2-23) 

I vAq ]  <  (y1+Y2+73)/(p3+P4), 

where  ( 1 , ( Y1+Y2 ) / (P3+P4 > > 1/P3 » 1/P4 )  =  (p^P2 ’P3 ) ‘  Here » 
we  can  choose  p^'s  to  obtain  the  most  convenient  forms  of 
the  right  hand  sides  of  (2.23)  for  the  later  uses.  Similar 
conditions  for  other  systems  are  given  in  Table  2.7. 

2.3.3  Sufficient  Conditions  for  Convergence  of  Y(t)  to  E 
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Table  2.7  Sufficient  Conditions  for  Convergence 


of  £(t)  to 

(2)  Y1+Y2n3  >0 

(3)  |2v2|A3j2+v3|A2|2|  <  (y1+Y2){(2y2+Y3)P3-2|vA0| }/(p3+p4), 

1 2vi !a3 I 2— V3 1 Ai I  2 1  <  (Y1+Y2){(2Y1+Y3)p4-2|vA0|}/(p3+p4), 
|vAq|  <  (y1+Y2+Y3)/(P3+P4), 

(Pi>P2»P3»p4) =  (l,(Y1+Y2)/(p3+P4),l/p3,l/p4) 

(4)  |v2|a3|2+v3|a2|2|  <y1(y2+y3), 

lvA0l<P2(Yl+Y2+Y3)’ 

(pi>p2,p3)  =  (1,P2|AQ| ,p3y1|a0| ) 

(5)  I A2 1  2 <  p3(Y1+Y3){(2Y2+Y3)/p3-!vA0| }/2, 

|vAq|  <  p3(y1+2y2), 

(ppP2.P3)  =  (P3<Y1+nr2)  ,i,p3) 

(6)  |v2|A3|2+v3|A2|2|  <min{p4(Y2+Y3)(Y1+Y2+Y3), 

Y  1{  (Y i+Y  2^  3 >  ( 1-P4  > * <P4/P2 )  I  vAQ  |  }  }  , 

| vAQ |  <  min { (p4/p2 ) ( Y1+Y2+Y3 ) , (Pj/p4 ) ( 1-P4 > ( Yj^Y^ > > > 
(p’,p^p’,Pi;)  =  (p4y1(y1+y2+y3),p2(y1+y2+y3),i,p4(y1+y2+y3)) 

(7)  I  I Ai I 2_vi lA2 1  2 1  <  (Y2+Y3){(y1+Y2+Y3)-(p1P3P4/(P1+P4))|vA0| }, 

I v3 1 A2 1 2_V2 ! A3 ! 2  I +(p2/p4)  lA2 1 2  <  IvAqI (Y2+Y3>/2, 

|vAq|  <min{(p3/2p1)(Y1+2Y2),(p3/(p4(p2+P3)))(Y1+Y3), 

(p1p4/<p3(p1+p4))Ky1+y2+y3)  , 

(P]_>P2 >P3>P45P3)  =  <PX  |vAQ  |  ,p2  |vAQ  1 2,p3  |vAq  |  ,p4  |vAQ  |  ,1) 
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In  this  section,  we  shall  show  that  if  £( t) -*■  Q ^  as 
t-*-°°,  then  the  corresponding  Y(t)  converges  to  £  as  t -*•<». 
It  is  sufficient  to  prove  the  following  statement: 


(S.2.2)  For  any  small  £  >  0,  there  is  a  6  >  0  such  that 


|£|  <  6  =>  inf  II Y-Y '  II  <  £  for  all  Y  e  {Y  :  |g|  <  6} 
Y'eE 

or 

d(£,£.)=sup  inf  ||  Y — Y f  ||  0  as  6 -*■  0 , 

0  YeEgY’eE 

where  E^  =  {Y  :  ||£  |  =  6}. 


The  above  statement  implies  that  if  £  is  close  to  0^,  the 
corresponding  set  E^  is  also  close  to  £  in  the  Y-space. 

Here,  we  give  a  proof  of  (S.2.2)  only  for  (2').  The 
proof  for  other  cases  can  be  established  in  a  similar  way. 
For  (2'), 

E6  =  {Y  :  C  =  y1y2y3  +  z1z2y3  +  z1y2z3+y1z2z3  =  6}. 

.  .  T 

Consider  a  line  parallel  to  (1, 0, 0 , 0, 0, G  )  and  through  a 

point  ( 0 , z-p y2 , z2 , y3, z3 )  .  This  line  intersects  E^  and  E  if 

/V  A  A  A 

y2y3  +  z2z3  i  0,  since  the  equations: 


,  A  A  AAA  A  a  A  .a  « 

yly2y3+ZlZ2y3  Zly2Z3+ylZ2Z3  =  6 » 
yly2y3  +  ZlZ2y3  +  Zly2Z3+ylZ2Z3  =  °*  ^ 


/V  A  /\  A 


has  a  unique  solution  (y^,y.p  if  y2y3  +  Z2Z3  ^  ^  ‘ 


(2.24) 


From  (2.24), 


,  -  .  .  A  A  A  A  , 

yl”yl  y2y3+z2z3  = 
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Hence , 


d(En  Sx)  <  S1  =  (Y  :  I  Y2y3+z2Z3  I  >  ^  ’ 

T 

Likewise,  by  using  a  line  parallel  to  ( 0 , 1 ,0 ,0 , 0 , 0 )  and 

^  ^  /v 

through  a  point  (y.^  »0  >y2  >z2  >Y3  ’z3 )  ,  we  have 

d(EnS2,E{oS2)<{!s,  S2  =  {Y  :  lz2y3+y3z2l  >6^}. 

Hence , 

d(Zn(S1uS2),  S6  n  (S1  uS2>)  <  s’5. 

Consider  the  sets: 


ig  =  tY  :  y2y3  +  z2z3  =  0,  z2y3+y2z3  = 

P2  4  {Y  :  |y2y3  +  z2z3|  =  S*5,  U2y3+y2Z3'  =  ^ ' 


A  line  parallel  to  (0,0,1,0,0,0)T  and  through  a  point 

(y^ z1,Q,z2,  y3,z3)T  intersects  F0  and  T2  if  y3*0,  since  the 

equation : 


y2y3±Z2y3  =  °’  ^3*  ^3=*^’ 


(2.25) 


has  a  unique  solution  for  (y2,y2).  From  (2.25), 


(y2-y^y3  = 


Hence , 

d(r0  n  r4,F2  n  r4)  <  616,  ru  r  (Y  :  |y3 \  >  Sk}. 

Since  c £,  we  have 

d(£  n  r 4 , Z ^  n  f4)  <  max{ 6^ , 6^} .  (2.26) 
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By  repeating  the  above  procedure  for  different  lines,  we 
obtain  results  similar  to  (2.26)  for  =  {Y  :  ) |  >6^}  and 
|  z  ^  |  >  5  ^ } ,  i  =  1,2,3.  To  sum  up, 

d(Z  n  (  u  (r4i  u  A4i) )  ,Z6  n  (  u(r4i  uA^)))  <  max{  6^ ,  6*}  . 
Since  Qy  e  E,  where  Qy  is  the  zero-vector  in  the  Y-space, 

d(z-(  u  (r4.  ua4.)),z6-(  u  (r4.  ua4.)))  <  bV. 

Hence , 

d(Z,Zg)  <  maxCe^A31* , 6^} . 

If  we  choose  a  6  such  that  max{6  6  ,6^}  <  e,  then  we  have  the 
desired  result.  Thus,  (S.2.2)  is  proved  for  (2'). 

As  mentioned  earlier  for  (2'),  5(t)  +  0  as  t  •+• «  for  any 
initial  conditions,  if  y1+y2+y3>0.  Hence,  by  (S.2.2),  if 
Y1+Y2+Y3>0)  Y(t)  converges  to  Z  as  t-*-°°.  For  other  systems, 
Y(t)  converges  to  Z  if  the  wave  intensities  satisfy  the 
conditions  in  Table  2.7.  We  do  not  know  whether  the  original 
systems  (3')-  (7')  satisfy  these  conditions,  that  is,  whether 
the  systems  have  attractors  satisfying  these  conditions.  On 
the  other  hand,  it  is  easy  to  find  natural  boundaries  for 
some  of  the  reduced  equations.  If  (vi»v2,v3)=  (1,1,1), 
(v'.v')  =  (1,1)  and  y^  >  0,  i  =  1,2,3  in  (3),  then 

I  s  -  dl  +  9Y3|vAQ|2,  (2.27) 

where  I  =  2x^ 2+x2 2+ ( x3+3vAq ) 2 ,  d  =  min{  2y^  >  2Y2  ,Y3  ^  ’  anc*  ^ 
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(2.28) 


(V1,V2,V3)  =  > 

i  <  -di  +  y3|vAq|2, 

where  I  = x^2+x22+2 (x3+vAQ ) 2 ,  d  =  minCy-^  ,y2  j  ?Y3  }  •  If 
(V1’V2 ,V3>  =  and  >  0»  i  =  1,2,3,  in  (4),  then 

I  <  -di  +  |  vAq2  J  2  /  Y-l  ,  (2.29) 

2  2  2  2 

where  I  =  2 (x1-(v^/y1)vAq  )  +x2  +x3  ,  d = min{ 2y1 ,Y2 >Y3 ) • 
Hence,  the  ellipsoid  on  which  I  =  0  is  a  natural  boundary. 
There  are  no  such  boundaries  for  other  cases.  But,  as  shown 
later  in  the  numerical  experiments,  some  pseudo-chaotic 
solutions  of  the  reduced  systems  are  bounded  in  certain 
regions .  Therefore ,  we  assume  that  the  reduced  systems  have 
attractors  in  which  the  trajectories  are  finally  trapped. 
Then,  a  sufficient  condition  for  Y(t)  to  converge  to  £  as 
t  -*•  °°  is  as  follows: 


(S.2.3)  Assume  that  a  subset  A  of  £  is  an  attractor 
of  equation  (2.10)  with  the  initial  condition  Y^e  £. 
Also  assume  that  A  has  a  smooth  boundary  3A  relative 

A 

to  £  and  that  at  any  point  Y^ e  3A,  F(Y^)  is  transverse 
to  3A.  Let  r  be  a  set  in  the  Y-space  such  that  if 
Y(t)  c  T  for  all  t  >  0 ,  then  £(t)  -*■  8^  as  t  00 .  Then,  if 
A  cr,  there  is  a  neighborhood  0  of  A  in  the  Y-space 
such  that  if  YQ  e  0,  y+(YQ)  c  0  and  Y(t)  =  /(t,Y0) 
converges  to  £  as  t-*-°°. 


The  sets  defined  in  (S.2.3)  are  sketched  in  Figure  2.3. 
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Figure  2.3  Sketches  of  Sets  in  (S.2.3) 


The  validity  of  statement  (S.2.3)  can  be  established  as 
follows.  Let  I  and  A  be  n-dimensional  surfaces.  Then,  3A 
is  a  (n-l)-dimensional  smooth  surface  by  assumption.  Let  us 
consider  a  set  Gyb  X  where  6  and  Gyb  an 

orthogonal  complement  of  a  subspace  which  is  a  translated 
tangent  plane  of  E  at  Yfa  to  the  origin.  Then,  Gy^  x  {Y^}  is 
( 6-n) -dimensional .  Since  3A  is  smooth,  there  is  a  neighbor¬ 
hood  N(3A)  of  3A  such  that  the  set 

NA  3A)  =  N  (  3  A )  n  u  Gv,  x  (Y,  } 

1  Yb  b 

Y,  e  3  A 
b 

is  smooth  (5-dimensional)  and 
WX(3A)  c  r. 

/v  OO  A 

Since  F  is  C  ,  there  is  a  neighborhood  W^OA)  of  3A  relative 
to  M^(3A)  such  that  F(Y)  is  transverse  to  N^(3A)  at  all 

A 

Y  €  NjOA).  Let  us  take  a  neighborhood  W£(E)  of  E  with  radius 

e  (i.e.,  d(3N  (E),E)=e)  such  that 
e 
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(N  (E)  n  N  OA) )  c  W2(3A)  . 

By  (S.2.1),  there  is  a  such  that 

d({Y  :  ZAY)  =  0 }  ,  { Y  :  $AY)  =  6})  <  e. 

Hence,  for  6<min{6^}, 
i 

Wfi<  E)  4  {Y  :  1  c  <  Y )  ||  ^  <  6}  c  Mg(E). 

Here,  N^CE)  is  a  neighborhood  of  I  and 
(Ng(E)  n  M^SA))  c  N  (3A)  . 

Hence,  y(t,YQ)  remains  in  T  as  long  as  y(t  ,YQ)  e  Ng(E) .  On 
the  other  hand,  Y(t,Yg)  remains  in  N^(E)  as  long  as  /(t,YQ) 
e  T,  since  ||  £  ( t )  ||  ^  decreases  if  Y(t)  Consider  an  open 

set  0  containing  A  such  that  30  consists  of  N,  (E )  n V^( 3A ) 
and  a  part  of  3Ng(Z).  Then,  O^r  and  if  YQeO,  then  V+(YQ) 
<=  0,  which  is  the  desired  result. 
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CHAPTER  3 


REDUCED  EQUATIONS 

In  this  chapter,  we  shall  choose  systems  from  Table 
2.4  which  may  have  periodic  or  chaotic  solutions  for  almost 
all  initial  conditions  in  the  X-space.  At  the  present  time, 
sufficient  conditions  for  ensuring  the  existence  of  chaotic 
solution  are  not  known  for  general  three-dimensional  ordinary 
differential  equations.  Many  numerical  experiments  showed 
that  the  existence  of  more  than  one  unstable  equilibrium 
points  could  lead  to  pseudo-chaotic  solutions.  Here,  we 
shall  restrict  our  attention  to  those  systems  whose  reduced 
equations  may  have  periodic  or  chaotic  solutions  globally 
(i.e.,  for  almost  all  initial  conditions)  in  the  X-space. 
Hence,  if  we  can  show  that  a  system  has  a  stable  equilibrium 
point  or  an  open  set  RQ  such  that  if  XQ  e  RQ,  JX(t)J  as 
t-*00,  where  X(t)  =  X(t,Xg),  then  the  system  is  discarded. 

In  Table  2.4,  and  v|  can  be  ±1.  The  presence  of 
both  positive  and  negative  (or  v|)  corresponds  to  the 
interaction  of  positive  and  negative-energy  waves.  It  is 
known  that  such  a  system  is  explosively  unstable,  that  is, 
one  or  more  of  the  wave  amplitudes  tend  to  infinity  in  finite 
time.  The  addition  of  linear  damping  terms  does  not  suppress 
such  a  fast  growth.  Therefore,  we  only  retain  those  cases 
in  which  all  the  v.'s  (and  v|'s)  have  the  same  sign.  Without 
loss  of  generality,  we  assume  that  all  the  v^'s  (and  v|'s) 


are  positive,  since  the  sign  of  each  (and  v|)  can  be 
altered  by  changing  the  signs  of  the  variables. 

The  procedure  for  eliminating  some  of  the  remaining 
cases  is  as  follows  :  (i)  We  assume  that  all  the  y^’s  are 
positive  ;  (ii)  If  there  are  more  than  one  equilibrium  points 
and  all  of  them  are  unstable  for  certain  value  of  parameters, 
then  the  case  is  retained  •,  (iii)  If  there  is  only  one  equi¬ 
librium  point,  or  if  at  least  one  of  the  equilibrium  points 
is  stable  for  any  values  of  parameters,  we  introduce  a  small 
linear  growth  term  to  the  system  by  changing  one  of  the  signs 

of  y . ' s  so  that 
1 

Y1  +  Y2  +  Y3  >  °  ’  (3.1) 

(iv)  If  the  condition  in  (ii)  holds  for  the  system  having 
one  linear  growth  term,  the  system  is  retained,  and  if  the 
condition  in  (iii)  holds,  the  system  is  discarded.  We  do 
not  introduce  more  than  one  negative  Yj's,  since  it  is  equiv¬ 
alent  to  an  equation  with  less  than  two  negative  Yj's  by 
the  substitution  t  -t  and  by  certain  variable  changes  ;  (v) 

If  there  exists  an  open  set  Rq  such  that  ||X(t)||  -► 00  as  t  00 
for  all  XQ  «  Rq*  where  X(t)  =  X(t,XQ),  then  the  case  is 
discarded . 

The  reduced  systems  discarded  by  the  foregoing  pro¬ 
cedure  may  have  local  chaotic  solutions.  The  original 
systems  corresponding  to  the  discarded  reduced  systems  may 
also  have  chaotic  solutions.  Moreover,  the  remaining  systems 
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may  not  have  chaotic  solutions.  For  the  class  of  three- 
dimensional  ordinary  differential  equations  such  that  the 
existence  of  only  one  equilibrium  point  implies  the  nonexist¬ 
ence  of  chaotic  solutions,  we  can  conclude  that  if  a  three- 
wave  interaction  system  has  global  periodic  or  chaotic 
solutions  in  the  X-space,  then  it  belongs  to  one  of  the 
remaining  cases. 

3.1  Analysis  of  Reduced  Equations 

3.1.1  System  (3) 

System  (3)  with  positive  y^'s  has  been  studied  by 
Pikovskii  et  al.  [12].  Here,  we  present  their  results  along 
with  more  detailed  analysis.  Without  loss  of  generality, 
we  can  assume  that  vAQ  <  0,  since  the  substitution  vAQ -► -vAQ 
is  equivalent  to  x^-*--x^,  i  =  l,3.  Then,  the  equilibrium 
points  are  PQ  =  {0^,}  and  =  (Xv>,  v  =  ±1: 

Xv  =  (vC(Y3/Yi)(  I  vA0 1  2_Y1Y2 I  vA0 1  ~ ' C I  vA0 1  2_yiy2 )?*} 

v[y1Y3(  |vA0[2-y1Y2)?5/{|vA0|-(  [vA0|2-Y1Y2^nis,  (  IvAq^-y^)^?  (3.2) 

if  they  exist.  At  PQ  and  P  ,  the  characteristic  equations 
of  dF/dX  for  (3)  are 

pQ(X)  =  (X+y3){X2+(y1+Y2)X+y1y2-|vA0|2}  =  0,  (3.3a) 

pv(x)  =  x3+(y1+y2+y3)x2+x[y2y3+y3y1|va0|/{|vA0|-(  lvA0i2_YiY2)l5} 

-(Y3/Yi)  ( I  vA0 1  2-Y1Y2  ^  I vAq  I  -(  I  vA0 1  2-YxY2 J*5}] 
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(3.3b) 


+4y3(|vAq|  -y1Y2)  =  0,  (3.3b) 

2  2 

respectively.  When  |vAQ|  <  Y-^Y2  =  [  vAq  I  >  Pg  is  asymptoti- 

2  ~  2 

cally  stable  and  P^  does  not  exist.  At  |vAq|  =  |vAg|  ,  P^ 

emerges  from  the  origin  and  (3.3b)  has  three  real  roots,  one 

2 

of  which  is  zero  and  the  others  are  negative.  When  AQ 

=  Aq  +  e,  where  £  is  positive  and  sufficiently  small,  equation 

(3.3b)  has  three  negative  real  roots. 

2  ~  2 

We  consider  the  case  where  A^  > A^  .  Here,  equation 
(3.3 a)  has  two  negative  and  one  positive  real  roots.  We 
assume  that  (3.3b)  has  one  real  root  -c  and  a  pair  of  complex 


roots  a±ib.  Then, 


=  c-2a’ 


^3^3 1  vA0 1 /{  I  vA0 1 "( I  vA0  *  -ylV  } 


-[ (y3/yx) ( I vAq I  2-Y1Y2 >**{  |vAq|-(|vA0| 2-Y1Y2 )ls} ]  =  a2+b2-2ac , 

4Y3( I vAo I 2_YiY2)  = c(a2+b2). 

2  ~  2 

Hence,  c  >  0.  We  assume  that  at  Aq  = Aq  ,  the  root  locus 
crosses  the  imaginary  axis.  Substituting  a=  0  into  (3.4) 
leads  to 


(3.4) 


( |  vAQ  |  2-Y1Y2  )1*{2 (y2+y3-y1)  |  vAQ  |  2+4y12Y2) 

=  IvA0I(y1(y1+y2)(y1+y2+y3)+2(y2+y3-y1)( |vA0 |2-y1y2>> • 


(3.5) 


Hence , 


2  2 

I vAq  |  <min{2Y1  Y2/(Y]_-Y2-Y3)  ,Y1(Y1+Y2)(Y1+Y2+Y3)/{2(Y1-Y2-Y3)}+Y1Y2} 


-  * 

'  i 


(3.6a) 


or 

I vAq | 2 >max{2Y1^Y2/(Y1-Y2-Y3)>Y1(Y1+Y2)(Y1+Y2+Y3)/{2(Y1-Y2-Y3)}+Y1Y2} 
£  |  vAg  |  2 .  (3.6b) 

Now,  we  consider  two  subcases  : 


Case  (i)  y^_  <  Y2+Y3  :  From  (3.5),  we  obtain  an  equation  for 

I vAo 1 2 : 

s(  |vA0|2)  4  4{(Y1-Y2)-Y32H  |vA0|2)2-[(y1-Y2)2(Y1+T2+Y3)2 
+8y1Y2{(y1-Y2)2+Y3(y1+Y2)}]  |vAq|2-16y13Y23  =0.  (3.7) 


Equation  (3.7)  has  a  positive  real  root  if  and  only  if  the 
first  term  is  positive  (i.e.,  >  Y2+Y3  or  Y^Y-^Yg). 

Since  Y-^  <  Y2+Y3  by  assumption,  equation  (3.7)  has  a  positive 

•  •  A  2  • 

real  root  if  and  only  if  Y2>Y1+Y3>  and  AQ  is  defined  by 

|vAq|2  =  [c1+{c12  +  64(y1Y2)3c2}15]/2c2,  (3.8) 


where 


2  2 

Cl=(Y'l~Y'2)  ^(Y1+T2+Y3)  +8y1Y2>+8y1Y2Y3(Y1+Y2)» 


c2  =  4{(y1-Y2)2-Y32}. 


Since  s(y1y2)<0,  yxY2  <  |vAQ|2. 

Case  (ii)  Y]>  Y2+Y3  :  Since  s ( | vAg | 2 )  <  0  and  s( | vAg | 2 )  >  0 , 

the  real  root  of  (3.6)  does  not  satisfy  (3.5).  Thus,  if 
Y2<Yi+Y3>  the  root  locus  does  not  cross  the  imaginary  axis. 
For  sufficiently  large  |vAg|,  equation  (3.3b)  has  one 
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positive  root  and  a  pair  of  complex  roots  since  the  coef¬ 
ficient  of  X  in  (3.3b)  is  positive  and  increases  monotoni- 
cally  as  |vAg|  increases.  From  (3.4),  we  obtain  an  equation 
for  2a: 

q(2a)  s  (2a)3+2(Y1+Y2+Y3)(2a)2+(Y1+Y2+Y3)(2a)+{2a+(Y1+Y2+Y3^ 

Cy2Y3+Y3Y1  I  vA0  1 7 {  I  vA0  I -( I  vA0  I  2-YlY2  )55}-[  ( y3/y1)(I  vA0  I  2~Y1Y2  ^ 
{|vA0|-(jvA0|2-Y1Y2^)^n-^Y3(|vA0|2-Y1Y2^  =  °>  (3.9) 

which  has  only  a  positive  (resp.,  negative)  real  root  for 
sufficiently  large  jvAQj  if  Y2  >  yi+y3  <resp.,  Y 2  <  Y^+Y 3 )  • 
From  the  above  results,  we  obtain  the  root  locus  of 
(3.3b)  as  shown  in  Figure  3.1.  The  origin  is  asymptotically 


Figure  3.1  Sketches  of  Root  Locus  of  (3) 


2  ~  2  .  . 

stable  if  Aq  <  Aq  ,  and  is  a  stable  node  m  a  biunstable 

plane.  If  y2  <yi+y3’  pv  asyniptotically  stable  for  any 
2 

value  of  Aq  .  If  y2  >Y1+Y3>  pv  1S  asymptotically  stable 
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when  | vAQ | 2  <  | vAQ | 2  <  | vAQ | 2 ,  and  when  | vAQ | 2  >  | vAQ | 2 ,  is 
an  unstable  focus  in  a  bistable  plane.  Thus,  system  (3) 
with  positive  Yj's  may  have  periodic  or  chaotic  solutions 
globally  if  Y2>yl+y3  and  AQ?  > ^Q2 ’ 

3.1.2  System  (4) 

Without  loss  of  generality,  we  assume  that  v>  0,  since 

the  substitution  v  -*•  -v  is  equivalent  to  x^  -►  -x^  and  x2  -*■  -x2  . 

Then,  the  equilibrium  points  are  Pn  =  {0V}  and  P  =  {X  }, 

u  a  vy  vy 

v,y  =  ±1: 

Xvy  =  (vy(y2/Y3)55{Y2(vpA02/(Y2/Y3>l5-Y1Y3)}15,  v^yA^/^/Yg^-Y^)*5, 

y(Y2/Y3)?5(vyA02/(Y2/Y3)i5-Y1Y3)/Y1)T>  (3.10) 

if  they  exist.  At  PQ  and  Pvy »  the  characteristic  equations 
of  dF/dX  for  (4)  are 

PQ(A)  =  (A+y1){A2+(y2+Y3)A+y2Y3-(vA(:|2/y1)2>  =  0,  (  3 . 11a) 

py(A)  i  A3+(y1+Y2+Y3H2+u(Y3/y2^(Y2/Y3+1)  I v|Aq2 

+4Y2(y|v|A02/(Y2/Y3)'5-Y1Y3)  =  0-  (  3  .  lib) 

We  shall  show  that  if  all  the  y^’s  are  positive,  then 
PQ  or  P  is  stable.  From  (3.11a),  if  |  v  |  AQ  2  <  y^  ( Y2  Y3  )**» 

PQ  is  stable.  If  |v|Aq  >yi^2y3^  *  there  exist  two  nontriv¬ 
ial  equilibrium  points  P^  and  P_1;]_  (i.e.,  y  =  l).  We  shall 
show  that  they  are  stable.  Since  |  v  |  AQ2  >  Y]_(y2Y3  )**,  p^a^>0, 
i=  0,1,2.  Hence,  if  (3.11b)  has  three  real  roots,  they  are 


negative.  Assume  that  (3.11b)  has  a  real  root  -c  and  a  pair 
of  complex  roots  a±ib.  Then, 


Yl+Y2+Y3  =  C"2a’ 

y(Y2+Y3)  |v|A02/(y2y3)1s  =  a2+b2-2ac,  • 

^^Y2Y3^I  v|Aq2-y1Y2Y3  =  c(a2+b2)/4. 

2  % 

Since  |v|Aq  >yi^y2y3'  and  y  =  l,  c  is  positive. 
Y1Y2Y3  =  ^Y2Y3/(Y2+Y3)-c/4}(a2+b2)-2acY2Y3/(Y2+Y3>, 
4Y1Y2/(Y1+Y2)  >c. 


(3.12) 


Furthermore 


Hence , 

c  =  Y1+Y2+Y3+2a<  4y1y2/(y1+Y2). 


This  implies  that  a<0.  Thus,  (3.11b)  has  no  unstable  com¬ 
plex  roots.  Therefore,  Pyl  is  asymptotically  stable  if  all 
the  Yj's  are  positive. 

If  Y,  <  0  or  y,  <  0,  P  does  not  exist.  This  case  is 
l  o  vy 

also  eliminated  in  Table  2.5. 


Now,  we  assume  that  y1<  0.  Equation  (3.11a)  has  one 

2  -v  2  a 

positive  and  two  negative  roots  if  AQ  <  AQ  %  (Y2Y3  I  Y-^/v  |  , 

and  has  two  positive  and  one  negative  roots  if  AQ  >  AQ  . 

2  ~  2 

Moreover,  if  A^  <  A^  ,  there  exist  four  nontrivial  equilib¬ 
rium  points  P  ,  v,y  =  ±1,  and  two  equilibrium  points  P^  and 

P-ll  lf  A0  >  A0  • 

We  first  consider  P  and  P_11  (i.e.,  y=l).  If 


11 


»•> 
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equation  (3.11b)  has  three  real  roots,  then  all  of  them  are 
negative  since  p^1  >0,  i =  0,1,2.  When  AQ  =  0  or  Aq  is 

sufficiently  large,  equation  (3.11b)  has  one  negative  real 
root  and  a  pair  of  complex  roots.  Moreover,  the  root  locus 

does  not  cross  the  origin,  since  p^CO)  t  0  for  any  values  of 

9  ... 

Aq  .  Suppose  that  (3.12)  holds.  Then,  c  is  positive,  since 

y  =  1  and  0.  We  consider  two  subcases  to  determine  when 

the  root  locus  of  (3.11b)  crosses  the  imaginary  axis: 

Cases  (i)  (Y2~Y3 ) 2- I Yx I (Y2+Y3 )  >  0  :  Let  a  =  0  in  (3-12)  >  then 

AC)2  =  4(Y2Y3)i5|Y1/v|  {(y2-Y3>2-|y1I  (Y2+Y3)}  =  Aq2  .  (3.13) 

2  ^2 

From  (3.12),  at  A^  = Aq  , 

da/dAQ2  =  |vj  {  Iy-^I  (Y2+Y3)-(Y2-Y3)2}/{2(y2y3)^(b2+c2)}  <  0,  (3.14) 

where 

£2  =  4y2y3Iy1|(y2+y3)/{(y2-y3)2-Iy1|(y2+y3)>, 
c  =  y2+y3-IyiI  • 

2  *  2 

Hence,  Hopf  bifurcation  occurs  at  Aq  =  Aq  . 

Case  (ii)  (Y2-Y3  )  2- |  Y-j_  I  ^Y2+Y3  )  <  0  :  Here,  (3.13)  does  not 

hold,  that  is,  the  root  locus  of  (3.11b)  does  not  cross  the 
imaginary  axis.  From  (3.12),  we  obtain  an  equation  for  2a: 

q(2a)  =  (2a)3+2(Y2+Y3-|Y1l )<2a)2+{(Y2+Y3-|Y1l )2 

+y(Y3/Y2)l5(Y2/Y3+D  |v|A02}+y(Y3/Y2)i5|v|A02{(Y2-Y3)2-|Y1|(Y2+Y3)}/Y3 

-t+|Y1lY2Y3  =  0,  (3.15) 
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where  y  =  1.  Equation  (3.15)  does  not  have  negative  real 

2  2 
roots  at  Ag  =0  and  for  sufficiently  large  Aq  .  This  means 

that  if  (3.11b)  has  a  pair  of  complex  roots,  their  real  parts 

.  2  .  2 
are  positive  at  =0  and  sufficiently  large  AQ  .  Thus, 

equation  (3.11b)  always  has  one  negative  real  root  and  a  pair 

of  complex  roots  whose  real  parts  are  positive. 

We  now  consider  ^  and  P  ^  ^  (i.e.,  y=  1),  which 
2  ~  2 

exist  only  if  Aq  <  Aq  .  Equation  (3.12)  with  y  =  -1  implies 

that  a  cannot  be  zero.  Hence,  the  root  locus  of  (3.11b)  does 

not  cross  the  imaginary  axis.  Equation  (3.12)  also  implies 

2 

that  the  real  parts  of  the  complex  roots  of  (3.11b)  at  AQ  =  0 
.  2  ~  2 

are  positive.  Moreover,  at  AQ  = AQ  ,  equation  (3.11b)  has 
three  real  roots  which  are  positive,  negative  and  zero. 

2  2  ~  2 

Hence,  b  =  0  in  (3.12)  for  a  certain  value  AQ  < AQ  .  From 
2 

(3.12),  Aq  is  uniquely  determined: 

A02  =  2(18  |Y1lY2Y3-a(Y2n3-|Y1l 

/C|v|[{6a-(Y2+Y3-|Y1|)KY2+Y3)+36Y2Y313,  (3.16) 

where  a  is  the  positive  root  of 

2(Y2+Y3>a3+{  (Y2+Y3^Y2+Y3- 1 Yx I  )+12Y2Y3>a2 

+8'Y2Y3(Y2+T3_lYll  )a_4lYlly2Y3(Y2+T3)  =  (3.17) 

Thus,  equation  (3.11b)  has  one  negative  root  and  a  pair  of 

2  2 

complex  roots  whose  real  parts  are  positive  if  Aq  <  Aq  ,  and 

2  2  ~  2 

has  one  negative  and  two  positive  real  roots  if  Aq  <  Aq  <  Aq  . 


From  the  above  results,  we  have  the  root  loci  shown  m 
Figure  3.2.  The  origin  P_,  is  a  stable  node  in  a  biunstable 


(Y2-Y3)  < 
lYl I (Y2+Y3) 


\  (Y2-Y3 >  > 

Pvl  ) Y 1  I C  Y  2  +Y  3 )  Pv.x 

Figure  3.2  Sketches  of  Root  Locus  of  (4) 


2  2 

plane  if  Ag  <  Ag  ,  and  an  unstable  node  in  a  bistable  plane 
.  2  2 

if  AQ  > Ag  .  In  case  (i),  Py^  is  an  unstable  focus  in  a 
.  2  ~  2 

bistable  plane  if  Ag  <  Ag  and  is  asymptotically  stable  if 
2  ~  2 

Ag  > Ag  .  In  case  (ii),  Py^  is  always  an  unstable  foci  in  a 
.  2  2 

bistable  plane.  If  An  < An  ,  P  .  is  an  unstable  focus  in  a 


2  2  ~  2 

bistable  plane  and  if  Ag  <  Ag  <  1  an  unstab3.e  node 

2^2  2 

in  a  bistable  plane.  Thus,  if  Ag  < Ag  ,  or  if  (y^-Yg)  < 

|YiKy2+Y3)  and  I  Y  Y  2 3  ((3.1)),  all  the  equilibrium 
points  are  unstable  and  therefore,  system  (4)  may  have  peri¬ 
odic  or  chaotic  solutions  globally. 

We  shall  show  that  system  (4)  has  an  attractor  contain- 
ing  P  .  Since  equation  (4)  is  invariant  to  the  interchange 
of  subscripts  2  and  3,  we  can  assume  that  y3 >  y2  without  loss 
of  generality.  We  assume  that  y3  >  Y2 ,  sance  Y3  =  Y2’  ^ 
becomes  a  two-dimensional  equation,  which  has  no  chaotic 
solutions.  Let  a(t)  =  x3 (t)/x2 (t)  ,  then 

a=  (y3-Y2)a-(x1-vAg2/y1) (1-a2) .  (3.18) 

At  x2  =  vx1,  v  =  ±1,  a  =  v  ( Y 3“Y2 )  •  Since  Y3>Y2»  then  the  sets 
V+  and  V_  defined  by 

V+ =  (X  :  x2  >  |x3| },  V_ =  {X  :  x2  < -|x3|},  (3.19) 

are  positive  invariant  sets  (i.e.,  if  X(x,XQ)e  V+  (or  V  ) 
for  some  x,  X( t ,XQ )  £V+  (or  V_)  for  all  t  >  x )  .  Let 
V  =  V+  u V_ .  Then,  it  is  enough  to  show  that  if  X(t,XQ)  «  -V 
for  all  t  e  [  0  ,°°)  ,  X(t)  =  X(t,Xg)  converges  to  V  as  t  ■+■ 00 .  Let 
8(t )  = x2 ( t ) 2-x3 ( t ) 2 ,  then 

B=2(y3x32-y2x22)=2(y3-y2)x32-y28.  (3.20) 

If  X  +  (Xg )  e  -V ,  then  | x2 ( t ) |  <  |  x3 ( t ) |  (or  6(t)<0)  and 
|  x  2  ( t )  |  >  (Y3/Y2)ls|x3(t)|  (or  8  ( t )  >  0)  for  t  >  0  .  Hence  , 
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that  is,  X(t)  =  X(t,Xg)  converges  to  V  as 


0(t)  0  as  t  -*■ 00 

t-*-°°.  Thus,  V  is  an  attractor.  Moreover,  ^  e  V+  and 

Pvl £  V-' 

3.1.3  System  (5) 

The  equilibrium  points  are  Pn  =  (0V)  and  P  =  {X  }, 

U  /v  V  V 

v  =  ±  1  : 

Xv  =  (y2>  v{-(Y2/Y3Hy1Y3+  I  vAq  |  2  )  (v2/y  3)vAQ)r ,  (3.21) 

if  they  exist.  At  Pq  and  P^,  the  characteristic  polynomials 
of  dF/dX  for  (5)  are 


pQ(X)  =  (X+2y2HX2+(y1+y3)X+y1y3+|vA0|2}  =  0,  (3.22a) 

pv(x) 4  a3+(y1+y3)x2+(i-2y2/y3)(y1y3+|va0|2) 

-2y2(y1Y3+|vAq|2)  =  0,  (3.22b) 


respectively.  We  shall  show  that  (5)  can  be  discarded. 

Equation  (3.22a)  implies  that  PQ  is  asymptotically 
stable  if  all  the  y^'s  are  positive. 

We  assume  that  y2  <  0.  If  (3.22b)  has  three  real  roots, 
they  are  negative  since  p^^(0)>0,  i=  0,1,2.  Assuming 
that  (3.22b)  has  one  real  root  -c  and  a  pair  of  complex 
roots  a±ib. 


Y-j+Yg  =  c-2a, 

(1-2y2/y3)(Y1Y3+|vA0|2)  =a2+b2-2ac,  - 

-2y2(y1Y3+|vAq|2) =c(a2+b2). 


(3.23) 
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Hence,  c  is  positive,  and  we  have  an  equation  for  2a: 

* 3 , n r/..  , _ .  ,2 . .  /. .  w..  . .  ■  i - 1 2, 


q(2a)  =  (2a)'3+2(Y1+Y3)(2a)^+{(Y1+Y3)2+(l-2-Y2/Y3)(Y1Y3+|vA0|2)}(2a) 


+(Yr2YlY2/Y3+Y3)(YlY3+lvA0l  5  =  °» 


(3.24) 


which  does  not  have  a  positive  real  root.  Therefore,  the 
real  parts  of  the  roots  of  (3.22b)  are  negative.  Thus, 
is  asymptotically  stable  if  y1<  0* 

Now,  we  assume  that  y^ <  0-  Equilibrium  point  exists 


only  if  |  vAq  |  <  )  y -^ !  Y 3  -  Equation  (5)  is  rewritten  as 

2 


•  > 

X1 

'-Yl 

O 

< 

> 

1 

X1 

X- 

vA_ 

X- 

3  J 

^  0 

-V3 

<  3 

x2 

=  - 

OJ 

>- 

N-/ 

-  x^ )x„ 

• 

-x. 


(3.25a) 


(3.25b) 


The  eigenvalues  of  the  square  matrix  in  (3.25a)  are  positive 

2 


and  negative  real  numbers  if  |vAQ|  <  Iy^Jy^  Hence,  there 
is  an  open  set  RQ  where  x1 <  0  and  |x2|  is  so  small  that  if 
*0  e  R0  and  =X(t,XQ),  |  x1(t)  |  ,  |x3(t)  |  -*•  ®  and  |  x2  ( t )  |  -*•  0 

as  t  -►  ® . 


Finally ,  we  assume  that  y3  <  0 .  If  | vAQ |  <  Y ^ I Y 3 i » 

2 

Pv  does  not  exist.  If  |  vA  Q  |  >Y1Iy3I  and  |  y  3  [  <Y-j_>  then 

2 

Py  is  asymptotically  stable.  If  |vAQ|  >  Y 1 1 Y  3  I  and  I  Y  3  I  >Yi; 
the  eigenvalues  of  the  square  matrix  of  (3.25a)  are  positive 
real  numbers.  Hence,  there  is  an  open  set  similar  to  Rg. 
Thus,  system  (5)  has  no  periodic  or  chaotic  solutions 
globally . 
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3.1.4  System  (7) 


The  equilibrium  points  are  Pn  =  {0V}  and  P  =  {X  }, 

U  A  V  V 


v  =  ±1: 

Xv=  ((y2Y3+|vA0|2)/y3,  v{-(Yi/y3)(y2y3+|va0|2)}\ 


(3.26) 


-v(vAq/y3)  {-(yx/Y3)  (V2Y3+ 1  vAq  | 2 )  }ls)T, 

if  they  exist.  At  Pg  and  P^ ,  the  characteristic  equations 
of  dF(X)/dX  of  (7)  are 


PQ(X)  ^  (X+y1)U2+(y2+Y3>X+y2Y3+|vA0|2}  =  0,  (3.27a) 

pv(x) £  x3+x2{(y1+y2+y3)-(y9Y3+I vAg | 2)/y3}+x{y1(y2+y3^ 

■3(yl/Y3)(Y2Y3+lvAo|2)}"2Yl(Y2Y3+lvAo|2)  =  °»  (3.27b) 


respectively . 

Equation  (3.27a)  implies  that  PQ  is  asymptotically 

stable  if  all  Y^'s  are  positive. 

We  assume  that  y2  <  0.  Equilibrium  point  Py  exists 
2 

only  if  | vAg [  <  |y2Iy3»  If  equation  (3.27b)  has  three  real 

roots,  they  are  negative  since  p  >0,  i  =  0,1,2.  Assuming 
that  (3.27b)  has  one  real  root  -c  and  a  pair  of  complex 
roots  a+ib, 


Yl+Y2+Y3“( Y2Y3+ I VA0 1 2 5 /Y3 = C‘2a ’ 

y1<y2+y3)-3y1Cy2Y3+  fvA0 1 2)/y3  =  a2+b2“2ac>  ' 

-2y1(y2Y3+|vA0|2)  = c(a2+b2). 


(3.28) 
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Hence,  c  is  positive,  and  we  have  an  equation  for  2a: 
q(2a) §  (2a)3+2{Y1+Y2+Y3-<Y2Y3+|vA0|2)/Y3}(2a)2 

+l’{Yl+Y2+Y3_(Y2Y3+lvA0^2)/Y3}2+Yl(Y2+Y3)_3Yl(Y2Y3+lvAo|2)/Y3](2a) 
+{Y1+Y2+Y3-^Y2Y3+IvAo|2)'/y3}{y1(y2+Y3)"3y1(y2Y3+!vAo!2)/y3} 
-2Y1(Y2Y3+!vA0|2)  =  0.  (3.29) 

If  |  vAq  |2<y3^Y3_2y2^/3,  q^^>0,  i  =  0,1,2.  Hence ,  the  real 

roots  of  (3.29)  are  negative,  and  the  real  parts  of  the 

2 

complex  roots  of  (3.27b)  are  negative.  If  |vAQj  > 

Y3 (y3~2y2 ) / 3 ,  a  cannot  be  zero  in  (3.28).  Moreover,  pv(0) 

2 

i  0  since  | vAQ |  <  |y2Iy3*  Hence,  the  root  locus  of  (3.27b) 

2 

does  not  cross  the  imaginary  axis.  Furthermore,  when  | vAq | 

=  y3 ( Y3“2y2 ) /3 ,  equation  (3.27b)  has  one  negative  real  root 

and  a  pair  of  complex  roots,  and  from  (3.29)',  the  real  parts 

of  the  complex  roots  are  negative.  Hence,  (3.27b)  has  three 

2 

stable  roots  for  | vAQ |  > y3 (y3-2y2  ) / 3 .  Thus,  is  asympto¬ 

tically  stable  if  y2  <  0. 

Then,  we  assume  that  y^  <  0 .  Rewriting  the  right  hand 
side  of  ( 7 ) , 


If  x^(0)  <0  and  |x^(0)j  is  sufficiently  large  and 


(3.30a) 


(3.30b) 


|x2(0)| 
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and  |  x  3  ( 0 )  |  are  sufficiently  small,  and  x2(t), 

x3  (t)  -*■  0  as  t  -*■  10 . 

2 

Now,  we  assume  that  y3<0.  If  |vAq|  <  Y  2  I  ^  3  I  ’  ^oes 
2 

not  exist.  If  |vAq|  >y2Iy3I  and  |  Y  3  I  <  Y  2  »  Pg  is  stat>le' 

2 

If  | vAq I  >  Y 2  I Y 3 |  and  | Y 3  I >  Y2  >  Pq  is  an  unstable  node  in  a 
bistable  plane.  Equation  (3.27b)  has  at  least  one  positive 
real  root.  Hence,  system  (7)  may  have  a  chaotic  solution, 

if  Y3  <  0,  |  Y 3  I  >  Y2  and  I  vAQ  j  2  >Y2Iy3I- 

3.1.5  System  (2)  and  (6) 

Since  system  (2)  can  be  obtained  from  (6)  by  the 
substitution  y^+avAQ  -*■  y^  >  we  consider  only  system  (6). 
Comparing  (2')  and  (6')  in  Table  2.3,  system  (6)  may  have  a 
chaotic  solution  only  if  Y-^+crvAg  <  0  and  y2»Y3  >  0.  Then, 
the  equilibrium  points  are  Pn  =  {0V>  and  Pv  =  (Xv  },  (k,u,v) 
e  K  =  {(1,1,1),(1,-1,-1),(-1,-1,1),(-1,1,-1)}: 

X<y  =  (k(y2Y3)15’  ^-(Y-j+ovAqJy^55,  v{-(y1+ovA0)y2}?5)T.  (3.31) 

The  characteristic  eqations  of  dF/dX  for  (6)  are 

pQ(X)  =  (X+y1+ovA0)(A+y2)(1+y3)  =  0,  (3. 32a) 

p(A)  ^  X3+X2(y2+Y3+Y1+ovA0)-4(y1+ovA0)y2Y3  =  0,  (3.32b) 

at  Pq  and  P^  ,  respectively.  Since  y1+ovAq<0,  equation 

(3.32b)  has  one  negetive  real  root  and  a  pair  of  complex 

roots  with  positive  real  parts.  Hence,  Pq  is  a  stable  node 

in  a  biunstable  plane  and  P^  is  an  unstable  focus  in  a 

k  y 
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bistable  plane. 

As  in  the  case  of  system  (4),  we  obtain  an  attractor 

V  =  V+  u  V  ,  where  V  and  €  V+  .  In  the  computer 

experiments,  we  did  not  observe  trajectories  such  that 

X  +  (Xg)€-V  and  X(t)=X(t,Xg)  converges  to  8V  as  t -*•“>.  We  do 

not  prove  the  nonexistence  of  such  trajectories,  but  we  can 

show  that  such  trajectories  are  highly  special  cases  which 

can  be  neglected  for  all  practical  purposes.  Since  (6)  is 

invariant  to  changing  the  signs  of  any  two  of  the  variables, 

it  is  enough  to  consider  the  set  {X  :  x3  >  0).  Let  us  define 

the  set  W  : 

tcp 

W  =  {X  :  Ikx-J  =  icx^,  |yx2|  =  yx2,  x3  >  0}  -  W,  K,y  e  {-1,0,1}} , 
where 

W=(  u  {Pv  })  u  {x_  ,x„,x--axis}. 

(k,V,v)«K  KV  1 

The  set  W  is  an  invariant  set  and  if  XQ*W,  X+(XQ)  n W =  $ . 

In  the  following,  the  trajectory  remaining  in  W  will  be 
neglected.  Consider  the  following  sets: 

S  4  {X  :  a  =  0}  n  (WQ0  -  V)  , 

W+(-)  =  {X  :  *  0  (a  <  0) }  n  (WQ0-V) , 

where  a(t)  =  x2(t)/x3(t)  and 

o=(Y3-y2)a  +  x^(l-a^).  (3.33) 

The  sketches  of  S  and  W+  are  shown  in  Figure  2.3,  In  W  , 
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X3 


Figure  3.3  Sketches  of  S  and  W+  in  (6) 

( <  » U )  =  (1,-1), (-1,1), 

jx1|=Kx1,  |x2!=<x2,  x2  <  0 .  (3.34) 

Hence,  the  trajectory  is  transverse  to  S  and  enters  W+  n  ^ 
In  W+  n  -1  <  a  <  0,  a  >  0 ,  x^>0  and  x^  >  0.  Hence,  from 

(3.33),  a  increases.  Thus,  the  trajectory  enters  WtnW^, 
Let  X(x1,XQ)  e  n  ( x1  ,x3 ) -plane  and  X(T1  +  e,XQ)  €  W+  n  W^, 
where  X(  t  ,XQ  )  e  W+  n  W.^  for  all  T  e  [t^  ,T^+e  ]  .  In  W+(iW^, 

0  <  a  <  1,  a  >  0  and  x^,  >  0 .  Hence ,  from  (3.33),  if  X(t,Xq) 
c  W+  n  W11  for  all  t  e  [T^+Gjtj, 

a(t)  >  (Y3“Y2)a(t)  >  ( Y 3+Y2 ( T i+e )  >  °*  (3.35) 

This  implies  that  the  trajectory  enters  W+ n  W  since  it 

is  transverse  to  the  ( x2 ,x3 )-plane .  Let  X(x2,Xg)e 
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WQ i n  (x2 ,x3)-plane  and  X(x2+s,Xq)  e  W+  n  W  where  X(x,Xq) 
e  W+nW  n  for  allxe[x2,x2  +  e).  In  W+  nW_^,  x^  <  0 ,  x^  <  0 , 
x2  >0,  x2  <  0,  x3  >  0  and  x^  <  0.  Hence,  if  X(x,Xq)  e  W+  n  W  ^ 
for  all  x  £  [x2+e ,t]  , 

x^t)  <  -(y1  +  avA0)x1(x2  +  e)  <  0, 

(3.36) 

x^(t)  < -y^x^(t) ,  i =  2,3. 

Consequently,  if  X(t,Xg)  e  W+  n  W  for  all  t  £  0,  then 
x^(t)  x2(t)  +  0  and  x3(t)-*-0  as  t -*•“>.  Otherwise,  the 

trajectory  enters  So  W^q.  Similarly,  if  X(x0,XQ)  e  S  nW_-^Q 
at  some  x^,  the  trajectory  converges  to  the  x^-axis  (where 
x^(t)-w)  as  t  -*■  «  or  enters  S  n  W^q  .  Thus  ,  if  X+  (XQ  )  c  WQ0  -  V  , 
then  as  t  •*■<*>,  the  trajectory  converges  to  the  x^-axis 
monotonically  after  some  or  no  oscillations  about  the  x3~ 
axis,  or  oscillates  about  the  x3~axis  for  all  t s  0 . 

Let  us  see  roughly  how  a  trajectory  oscillating  about 
the  x3~axis  for  all  t>0  behaves.  Rewriting  (3.20)  for  (6), 

6 =  2 (y3x32-y2x22 )  =  2 ( Y3-Y2 )x32-2y2S •  (3.37) 

Integrating  (3.37), 

0  <  /QX3(x)2dx  =  (B(t)-B(0)  +  2Y2/QS(x)dx)/2(Y3-Y2^ 

<  -B(0)/2(y3-Y2) ,  (3.38) 

for  all  t  since  B(t)  <  0  for  all  t.  By  studying  F(X)  of  (6) 
carefully,  we  know  that  x2(t)  takes  a  larger  value  than 


{-(y^+ovAq )Yj}  at  each  time  when  the  trajectory  encircles 

the  x^-axis.  Therefore,  for  any  small  e  >  0 ,  the  time 

2 

duration  for  which  Xg(t)  >e  converges  to  zero  as  t  -*■». 

2 

Thus,  x^(t)  behaves  as  a  train  of  pulses  whose  heights  are 
larger  than  {- ( y^+cvAq ) Y 3 and  whose  widths  converge  to 
zero  as  t  -*■ 00 .  The  above  arguments  also  hold  for  x^. 

A  trajectory,  which  lies  in  W+  n  W  ^  or  W  n  ^1  1  ^or 
all  t  >  for  some  and  converges  to  the  x^-axis  monotoni- 
cally,  is  unstable  in  the  sense  that  a  small  perturbation  of 
a  can  shift  X(t,Xg)  into  V  or  cause  X(t,XQ)  to  oscillate 
about  the  Xg-axis.  For,  a  =  p(Y3-Y2^  on  9V  (|y|  =  yy)  and  the 
trajectry  is  transverse  to  S,  and  3  V  and  S  approach  to  each 
other  as  |x^|  A  trajectory  which  oscillates  about  the 

x3-axis  for  all  t  ?  0  is  also  unstable  in  the  sense  that 
X(t,XQ)  is  shifted  into  V  by  small  perturbation  of  a.  For, 
|x^(t)|  is  larger  than  {  -  ( Y^+crvAg  )  y^  ,  i=  2,3,  at  each  turn 
around  the  Xg-axis,  and  X(t)  =X(t,XQ)  converges  to  3V  as 
t Hence,  the  trajectories  which  remain  in  -(W  uV)  for 
all  t^  0  can  be  neglected  for  all  practical  purposes.  In 
fact,  such  trajectories  were  not  observed  in  computer 
experiments . 

We  note  that  (6)  is  the  only  system  which  explicitly 
includes  a  =  sin( 0 ^ ( 0 ) +9 ^ ( 0 ) ) .  If  |vAq|<|y^|,  P^y  exists  and 
is  unstable  when  y^  <  0 ,  and  Pq  is  stable  when  y ^  >  0  •  If 
| vAg |  >  Iy-J  ,  P^g  exists  and  unstable  when  sgn(avAg)  =  -l,  and 
Pg  is  stable  when  sgn(ovAg)  =  1.  We  assume  that  vAg  >  I Y ^ I 
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and  that  the  trajectories  converge  to  Z  as  t  -*■  <*> .  The 
question  is  which  of  the  reduced  equations  with  o  =  1  or  -1 
describes  the  asymptotic  behavior  of  (6')  in  the  limit  t  ->•<» . 
Obviously,  the  equation  with  a  =  1  is  not  the  one,  since  the 
origin  of  the  Y-space  is  unstable  if  vA^  >  |  -y  I  *  ^  is 

possible  that  either  one  of  them  does  not  describe  the 
asymptotic  behavior  of  (6')  in  the  limit  t-*-00,  but  the 
trajectory  of  the  original  system  approaches  in  turn  either 
of  the  subsets  of  1  corresponding  to  the  X-spaces  of  the 
systems  with  a  =  1  and  -1. 

3.2  Physical  Interpretation  for  Reduced  Systems 

In  this  section,  we  shall  try  to  give  some  physical 
interpretation  for  the  results  presented  in  the  previous 
section.  The  equations  for  frequency  matching  and  their 
corresponding  decay  diagrams  are  shown  in  Table  3.1. 

3.2.1  System  (2) 

In  the  decay  process  of  (2),  wave  1  decays  into  waves 

2  2  2 

2  and  3.  Let  us  define  I =  2x^  +x2  +x3  .  For  convenience, 
we  call  I  the  energy  of  (2),  though  it  is  not  exactly  the 
wave  energy  of  (2).  Taking  the  time  derivative  of  I  gives: 

I  =  -  2(y1x12 + y2x22  + y3x32) .  (3.39) 

We  assume  that  Y-j.  <  0.  The  constant  I  and  1=0  surfaces  are 
sketched  in  Figure  3.4,  where  they  correspond  to  the 
ellipsoid  E  and  the  elliptic  cone  C,  respectively.  If  X(t) 
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Table  3.1  Frequency  Matching  and  Decay  Diagrams 


X3 


Figure  3.4  Dynamics  of  Energy  of  (2) 

is  inside  (resp.,  outside)  C,  I  increases  (resp.,  decreases). 
That  is,  when  the  amplitude  of  wave  1  is  large,  the  linear 
growth  of  wave  1  is  dominant  so  that  the  energy  of  the  whole 
system  increases.  When  the  amplitude  of  wave  2  or  3  is 
large,  the  linear  damping  of  2  or  3  is  dominant  so  that  the 
energy  of  the  whole  system  decreases.  Thus,  the  energy  of 
the  waves  are  not  balanced  in  a  static  manner,  but  may 
oscillate  in  a  chaotic  manner. 

As  mentioned  in  Section  2.1,  if  y2  <  0  (°r  Y3  >  0), 
there  is  no  chaotic  solution.  For,  if  y?  <  0, 

(d/dt)(|A1|2+lA3|2)  =  -y1!a1|2-y3|a3|2, 

(3.40) 

(d/ dt) ( | Ax | 2- | A2 | 2 )  =  -y1|A1|2-|y2I  | A2 I  2 ,  . 

and  therefore  ,  |  A^(  t )  |  and  [A3(t)|-*-0and  |  Aj  ( t )  |  •*>  • 

as  t -*■<*>.  This  phenomenon  does  not  depend  on  the  magnitude 

of  Y^'s,  since  they  are  normalized  by  the  substitution 
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| 1 1  -*■  t ,  i  s  1,2  or  3.  From  the  physical  view  point,  these 
results  seem  to  be  peculiar.  For,  system  (2)  describes  a 
three-wave  interaction  including  both  decay  and  fusion 
processes,  and  it  could  be  expected  that  the  energy  of  the 
linearly  growing  and  damped  waves  are  balanced  as  in  the 
case  where  y^  <  0.  The  above  results  show  that  when  the 
decaying  wave  1  is  linearly  growing  (y^<0),  its  energy  is 
transferred  to  the  damped  waves  2  and  3.  On  the  other  hand, 

when  wave  2  is  linearly  growing  (.y ^  <  0),  its  energy  is  not 

transferred  to  other  waves.  It  is  suggested  that  when  the 

linearly  damping  or  growth  terms  are  introduced,  even  if  the 
damping  or  growth  rates  are  very  small,  system  (2)  describes 
only  the  irreversible  decay  process.  It  should  be  noted 
that  the  nonlinear  coupling  terms  in  (2)  does  not  possess 
any  mechanisms  for  such  an  irreversible  prpcess. 

3.2.2  System  (3) 

As  mentioned  earlier,  in  the  decay  process  of  (3),  an 
external  wave  with  constant  amplitude  decays  into  waves  1 
and  2  which  have  another  resonant  wave  3.  When  the  external 
wave  is  weak  (|vAQ|  "t*le  ^inear  damping  of  three 

waves  is  dominant.  Therefore,  the  energy  supplied  by  the 
external  wave  is  dissipated  through  the  linearly  damped 
waves.  If  the  external  wave  is  strong  (  |  vAq  |  the 

input  energy  is  supplied  at  a  higher  rate  than  the  dissipa¬ 
tion  rate  when  the  amplitudes  of  damped  waves  are  small. 
External  energy  and  energy  dissipation  are  finally  balanced 
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when  the  amplitudes  of  waves  are  large.  When  wave  2  is 

highly  damped  linearly  (y2>y^+Y3),  the  external  energy 

supply  and  energy  dissipation  are  balanced  in  a  static 

manner  regardless  of  the  magnitude  of  the  external  energy. 

That  is,  the  energy  of  each  wave  reaches  a  constant  value. 

When  wave  2  is  damped  slowly  (y2  < y^+Yg)  and  the  external 

2  ~  2 

wave  is  not  strong  (A^  <  Aq  ) ,  the  wave  energy  is  also 

balanced  in  a  static  manner.  When  the  external  wave  is 
2  *2 

strong  (Aq  >  Aq  ) ,  the  energy  of  the  waves  may  oscillate  in 
a  chaotic  manner. 

3.2.3  System  (4) 

In  the  decay  process  of  (4),  the  external  wave 
generates  a  harmonic  wave  1  which  decays  into  waves  2  and  3. 
If  all  the  waves  are  linearly  damped  and  the  external  wave 
is  weak  ( j v ]AQ  < y^(y2Y3 )^) ,  the  energy  is  dissipated  more 
rapidly  than  supplied.  When  the  external  wave  is  strong 
( |v|AgZ  > y1(y2Y3)^) ,  the  energy  supply  is  at  a  higher 
(resp.,  lower)  rate  than  the  dissipation  if  the  wave  ampli¬ 
tudes  are  small  (resp.,  large).  Thus,  the  wave  energy  is 
balanced  in  a  static  manner  at  certain  values  of  amplitudes. 
If  wave  2  (or  3)  is  linearly  growing,  wave  2  (or  3)  grows  to 
infinity  as  t  -*■  <*>  regardless  of  the  magnitudes  of  the  growth 
or  damping  rates.  This  phenomenon  is  similar  to  that  of  (2) 
discussed  in  Section  3.2.1. 

If  wave  1  linearly  growing  and  if  there  is  no  external 
wave,  the  system  is  identical  to  (2).  When  wave  1  grows 
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2 

rapidly  C  I  Y I  ^ “Y 2 3 ^  >  '  or  t*le  external  wave  is  weak 

2  ~  2 

(Aq  < Aq  ),  the  system  behaves  in  a  similar  manner  as  (2). 

2 

If  wave  1  grows  slowly  (|y1|Cy2+73)  >  ( ^ 2~^ 3 ^  ^  an^  t^ie 
•  2  ^  2 

external  wave  is  strong  (Aq  >  Aq  )  ,  the  energy  supply  is 
dominated  by  the  external  wave  and  the  wave  energy  is 
balanced  in  a  static  manner. 

3.2.4  System  (5) 

As  shown  in  Table  3.1,  we  have  two  cases.  In  the 
decay  process  of  the  first  case,  the  external  wave  generates 
wave  3  by  interacting  with  1  which  is  a  harmonic  wave  of 
wave  2.  In  the  decay  process  of  the  second  case,  the 
external  wave  generates  wave  1  by  interacting  with  wave  3, 
and  wave  1  generates  a  subharmonic  wave  2.  If  all  the  waves 
are  linearly  damped,  the  energy  of  the  whole  system  is 
dissipated  regardless  of  the  external  wave.  If  wave  2  is 
linearly  growing,  the  wave  energy  is  balanced  in  a  static 
manner  regardless  of  the  amplitude  of  the  external  wave. 

If  wave  3  is  slowly  growing  linearly  C  ( y 3  [  <y2^  an<^  ^e 
external  wave  is  weak  (|vAq|  <  I T 3  I  5 »  wave  1  and  3  grows 
to  infinity  and  wave  2  dies  away.  If  the  external  wave  is 
strong  ( | vAq |  >  Y ^  I Y 3  I ^ »  wave  3  is  strongly  coupled  with 

linearly  damped  waves  1  and  2  and  consequently,  all  waves 
die  away.  If  wave  1  is  growing  rapidly  ( | y  3  I  <Y^^j  wave  1 
and  3  grow  to  infinity  and  wave  2  dies  away  regardless  of 
the  amplitude  of  the  external  wave.  If  wave  1  is  growing 
slowly  (  |  y-l  {  <Yg)  and  the  external  wave  is  strong 
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2 

(  |  vAq  |  >  I  Y-jJ  Y3  ) »  waves  die  away.  If  the  external  wave 

2 

is  weak  (|vAg|  <  J  y I Y 3 )  »  the  wave  energy  is  balanced  in  a 
static  manner,  or  wave  1  and  3  grows  to  infinity  and  wave  2 
dies  away.  The  mode  of  the  actual  behavior  depends  on  the 
initial  conditions  if  wave  2  is  highly  damped  (2y2>Y3)> 
and  if  2y2  <  Y3»  it  depends  on  the  y^'s  and  vAq  in  a  more 
complicated  way. 

3.2.5  System  (6) 

In  the  decay  process  of  (6),  the  external  wave 
generates  a  subharmonic  wave  1  which  decays  into  waves  2 
and  3.  We  can  make  the  same  argument  as  we  used  for  (2)  by 
replacing  y^  with  y-^+ovAq.  If  y^  >  °>  and  the  external  wave 
is  weak  (|vAQ|  <y^)t  all  waves  die  away.  If  y^  >  0  and  the 
external  wave  is  strong  ( [ vAQ |  >Y1)*  there  may  be  a  chaotic 
solution.  If  y-^  <  0>  there  may  be  a  chaotic  solution  regard¬ 
less  of  the  external  wave. 

3.2.6  System  (7) 

Here,  we  have  two  distinct  cases.  In  the  decay 
process  of  the  first  case,  the  external  wave  generates  wave 
2  by  interacting  with  wave  3  and  wave  2  generates  a  harmonic 
wave  1.  In  the  decay  process  of  the  second  case,  the 
external  wave  generates  wave  3  by  interacting  with  wave  2 
which  is  a  subharmonic  wave  of  wave  1.  If  all  the  waves  are 
linearly  damped,  the  energy  of  the  system  is  dissipated 
regardless  of  the  external  wave.  If  wave  2  is  linearly 
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growing  and  the  external  wave  is  weak  ( | vAq |  <  min{ | y2 I Y3 > 

YsCYi+Ys)}),  the  wave  energy  is  balanced  in  a  static  manner. 
If  the  external  wave  is  strong  ( | vAq j  >  j  y 2  !  Y  3 ) 9  all  waves 
die  away  since  wave  2  is  coupled  well  with  the  linearly 
damped  wave  3  through  the  external  wave.  If  wave  1  is 
linearly  growing,  wave  1  tends  to  infinity  and  wave  2  and  3 
die  away  regardless  of  the  external  wave  amplitude.  If  wave 
3  grows  rapidly  C  |  y 3  |  >  Y2>  an^  "the  external  wave  is  strong 
( I vAQ |  > y2Iy3I  ) >  the  wave  energy  may  oscillate  in  a  chaotic 
manner. 

3.3  Systems  with  Parametric  Instability 

Now,  we  have  six  reduced  systems  which  have  periodic 
or  chaotic  solutions  for  almost  all  initial  conditions. 

These  systems  are  listed  in  Table  3.2,  where  all  the  y^'s 
are  positive  except  for  (6).  The  parameters  of  each  system 


satisfy 

the  following  condition: 

(3) 

Y2  >  Yi 

+Y3,  A02  > A02  ((3.8)), 

(3.41) 

(4) 

2 

V<A 

(y2-y3 

Q2  ((3.13)), 

2 

)  <Yj(y2+Y3)>  "^l<"^2+^3’ 

(3.42) 

(6) 

Yx  >  0  > 

ctvAq  <  -y x  j  or 

(3.43a) 

Yx  <  0, 

ovAq  >  0 , 

(3.43b) 

(7) 

CM 

O 

1 

>y2y3,  y3>y2. 

(3.44) 
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Table  3.2  Systems  for  Chaotic  Solutions 


(1* ) 

Al"YlAl  =  -i{A22+AiCviiilAi  1 
A2+Y2a2  =  -i{A1A*+A2(V2iitA: 

2+V122 Ia2 |2)}> 

|2+V222IA2|2)}’ 

(1) 

Xi  =  y1x1-6x2+2x1x2-x2( (V1U 
2 

x2  =  ‘'rix2+<Sxi+x3~2xi  +X1^V2 
x3  = -2y2X3-^x2X3 ’  x3  -  0 

-2V2n)(x12+x22)+(V122-2V222)x3}> 
.11~2V211)(x1  +x2  )+(V122-2V222)x3^’ 

(2') 

A1  =  Y1A1-iA2A3, 

A2  =  -YjAj-iAjA*, 

A3  =  -Y3A3“iAlA2 

(2) 

x1  =  Y1x1-x2x3, 

x2  =  -Y2x2+x1x3, 

x3  =  ~Y  3X3+X1X2 

(3') 

A1  =  -Y1A1-i(A2A3+vA(JA*) , 

A2  =  “Y  ^2  “  i ( AxA*  +vA0Al  ^  ’ 

A3  =  ’Y3A3-iAlA2 

(3) 

X1 = “Y1X1-X2X3-VA0X2 ’ 

X2  =  -Y2X2+X1X3_VA0X1’ 

X3  =  _Y3X3+X1X2 

(4’) 

A1  =  YlAl"i(A2A3+vA02 ) ’ 

A2  =  _Y2A2_iAlA3’ 

A3  =  -Y3A3-iAlA2 

(4) 

x2  =  -Y2x2-x1x3+(vA02/Yi)x3, 

2 

x3 = -YgXg-x^Xj+CvAQ  /Yj)x2 

(6’) 

A1  =  ~YlAl~i(A2A3+vA0Al)  ’ 

A2  =  *Y2A2-;LA1A3’ 

A3  =  -Y3A3-iAlA2 

(6) 

Xl=  (-Y1-ovA0)x1-x2x3, 

x2  Y2X2+X1X3’ 

X3  =  -Y3X3+X1X2 

(7* ) 

Al  =  _YlAriA2  ’ 

A2  =  ■Y2A2"i(AlA2+vA0A3>’ 

A3  =  Y3A3~ivA0A2 

(7) 

%  =  -YlVX22’ 

x2  =  ~Y2x2+x^x2+vAqX3 , 

X3  =  y3x3"vA0X2 
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The  systems  given  in  Table  3.2  except  for  (7)  have  an 
unstable  focus  in  a  bistable  plane.  All  three-dimensional 


ordinary  differential  equations  having  pseudo-chaotic 
solutions  so  far  have  this  type  of  equilibrium  point. 
Actually,  we  did  not  observe  any  pseudo-chaotic  solutions 
of  (7). 

Systems  (1),  (2),  (4)  and  (7)  have  linear  growth  terms 
which  cause  the  equilibrium  points  to  be  unstable.  Such 
models  are  not  widely  accepted,  since  the  linear  growth 
terms  cannot  be  readily  derived  from  the  fundamental  equa¬ 
tions  describing  wave-wave  interactions  in  plasmas.  On  the 
other  hand,  introducing  linear  damping  terms  is  natural 
since  any  collective  motion  of  plasmas  is  damped  by  particle 
collisions.  System  (3)  has  only  damping  terms,  and  so  does 
system  (6)  if  y^  >  0  ((3.43a)).  Hence,  in  what  follows,  we 
consider  only  these  two  cases  where  the  instability  is  due 
to  an  external  wave  (parametric  instability). 

Now,  we  combine  the  sufficient  conditions  for  the 
convergence  of  the  trajectories  to  Z  (Table  2.7  and  (S.2.3)) 
with  the  conditions  for  the  existence  of  global  periodic  or 
chaotic  solutions  in  the  X-space  ((3.41),  (3.43a)). 

3.3.1  Sufficient  Conditions  for  Convergence  of  §(t)  to  0^ 

We  now  recall  the  sufficient  conditions  for  £(t)  to 
converges  to  0^  as  t  ■*■<*>.  Rewriting  the  equation  (3)  in 
Table  2.7, 
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4 


lA2 1 2+2 lA3 1 2  <  (y1+y2){(2y2+Y3)p3"2IvAoI >/(P3+P4)» 

I  I  a2  I  2-2  iA3 1 2  I  <  (Y1+Y2H(2Y1+Y3)p4-2|vA0|}/(p3+p4),  • 

|vA0|  <  (y1+Y2+Y3)/(P3+P4)- 

For  system  (6),  by  choosing  p2  and  p4  satisfying 
p4/p2  =  |vAq|/(y1+y2+y3), 


(3.45) 


the  right  hand  side  of  (6)  in  Table  2.7  is  maximized: 

Ia2  1 2+  lA3 1 2  <  Yi(Y2+Y3H1-|vAo|2/(Yi+Y2+y3)2}'  (3.46) 

We  denote  the  set  whose  elements  satisfy  (3.45)  (or  (3.46)) 
by  F.  Equation  (3.45)  and  (3.46)  implies  that  the  volume  of 
T  decreases  as  |vAQ|  increases  to  its  maximum  value  if  all 
the  y^’s  are  fixed.  We  also  know  that  as  y^>  i  =  1,2,3, 
becomes  large,  the  maximum  value  of  |vAQ|  and  the  volume  of 

/v 

T  increase.  As  y^  "*■ 00 ,  i=l,2  or  3,  and  max(Y^  }  /  |  vA^  |  ■+  ® , 

T  is  extended  to  the  whole  Y-space. 

3.3.2  Attractor  in  X-space 

We  obtain  an  attractor  of  (3)  smaller  than  that 
described  by  (2.27)  in  a  manner  similar  to  that  given  by 
Treve  [27]  for  the  Lorenz  system.  Equation  (2.27)  was 
obtained  as  follows: 


(d/dt){2x12+x22+(x3+3vAQ)2} =  -4Y1x12-2Y2x22-Y3(2x32+6vAC)x3) 

<  -2y1(2x12)-2y2x22-Y3(x3+3vAo)2+9y3IvAo|2 
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< -d{ 2x12+x2 2+(x3+3vAq ) 2  >+9y3 J  vAq j 2 ,  d  § min{ 2y1 , 2y2 ,y3 } . 


(3.47) 


Hence,  the  set  fig: 

rz0  ^  {X  :  2x12  +  x22  +  (x3  +  3vAq)2  <  9  (y  3 /d)  |  vAq  |  2  }  (3.48) 

is  an  attractor.  Suppose  that  y3  =  d,  and  let  y^  and  y2 
increase  while  j  vAg  |  is  fixed.  Then,  from  (3.45)  and  (3.48), 
fig c  T  is  satisfied  for  sufficiently  large  y^  and  y^.  Hence, 
by  (S.2.3),  there  exists  a  set  0.  Likewise,  for  suffi¬ 
ciently  large  y3  and  fixed  y3>Y2  and  |vAg|,  or  for  suffi¬ 
ciently  small  |  vAg  |  and  fixed  y^,y2  and  Y3>  there  exists  a 
set  0 . 

A  2  2 

Taking  the  time  derivative  of  1^ =  x^  +(x3+vAQ)  and 
I2  =  x22-(x3-vAq)2 , 

■^1  -  "dlIl+73  (vA0  I  2’  d1  =  min{2y1,y3},  (3.49a) 

~2Y2X22+Y2(x3'vA0)2_Y3lvAo|2’  if  2y2  >  Y3 ’ ' 

I2  >  •  ■  (3.49b) 

•  — Y 3 1 2  — Y 3  I vA0 I  2 ’  if  2Y 2  “  Y 3 • 

Thus,  the  sets  fi^  and  fi2  : 

fi-j^  §  {X  :  x]_2  +  (x3+vAq)  2  (Yg/dj^)  |  vAQ  |  2},  (3.50a) 

n2  §  {X  :  -x22+(x3-vAg)2  <  |vAg|2},  (3.50b) 

are  attractors  (it  was  assumed  that  2y2  < Y3  for  fi2).  Hence, 

^  a  ~  * 

A  =  fig  n  is  an  attractor,  and  A  =  fig  n fi^  n  ^2  is  an  attractor 
if  2y2  <  Y3 •  Figure  3.5  gives  sketches  of  A  and  A.  Here, 
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X 


% 

> 


Figure  3.5  Attractors  A  and  A  of  (3) 

if  2y2  5  Y3>  the  dint  at  the  bottom  appears  when  C y 3 / d )  >  9. 

3.3.3  Phase  Locking  and  Chaotic  Solutions 

Here,  we  combine  equation  (3.45)  with  Figure  3.5  to 
see  when  A c  f  is  satisfied.  If  Acf,  then  by  (S.2.3),  there 

A  A 

exists  a  set  0.  If  A c  r  is  not  satisfied,  the  sufficient 

A  A 

condition  A  c  T  for  the  existence  of  0  is  too  strong  or  there 
exists  no  0, 

Suppose  that  all  the  Y^'s  are  fixed.  As  mentioned  in 
the  previous  section,  A  c  f  is  satisfied  and  there  exists  a 
set  0  for  sufficiently  small  |vAq|.  In  this  case,  phase 
locking  is  of  no  importance,  since  the  trajectory  converges 
to 


We  consider  under  what  conditions  A  c  T  is  net 
satisfied.  From  (3.45), 


A2|2+2|A3|2  <  (2y2+Y3)(y1+y2), 

(3.51a) 

|A1 | 2-2 |A3 | 2 |  <  (2y1+y3)(y1+Y2)  , 

(3.51b) 

vAq|2  <  (2min{y1,Y2}+Y3) (y1+Y2+Y3)  . 

(3.51c) 

A  /V 

Hence,  a  necessary  condition  for  4  =  F  is  as  shown  in  Figure 
3.6.  Here,  the  projection  of  A  onto  the  ( x2 , x3 ) -plane  and 


Figure  3.6  Necessary  Conditions  for  Acf 


the  (x^ ,x3 )-plane  satisfy  (3.51a)  and  (3.51b),  respectively. 
Then,  from  (3.48)  and  (3.51), 

|  vAq  |  2  <  (2y2+y3)  (y^+y2)/2  (9(Y3/d)  +  (y3/d1)}/2-{l-4(y3/d)15}  ,  (3.52a) 


I vAn | 2 <  (2y. +y.) (y. +y7 )/{!+( y,/d1 fi}2 . 


(3.52b) 


Since  d<d^,  from  (3.52), 

| vAg |^<(2y2+y3)(y1+Y2)/16‘  (3.53) 

As  mentioned  in  Section  3.3.1,  equation  (3)  has  two 

2  ~  2 

nontrivial  equilibrium  points  P^'s  if  Ag  >  Aq  .  If 
y2  >  y1+y3’ 

|vA0|2  =  Y-,^2  >  (2Y2+Y3)  (y1+Y2)/16  •  (3.54) 

2  ~  2 

Hence,  if  Y2>Y1+Y3’  ^3-52^  does  not  hold  for  AQ  > AQ  . 

In  other  words,  when  system  (3)  has  periodic  or  chaotic 

2  ~  2  a  ^ 

solutions  globally  (AQ  >  AQ  ) ,  the  sufficient  condition  A  c  r 
is  too  strong  so  that  we  cannot  determine  whether  there 
exists  a  set  0. 

It  is  likely  that  there  may  exist  an  attractor  smaller 

A 

than  A,  since  equation  (3.52)  is  determined  at  points  P  and 
Q  in  Figure  3.6.  Equation  (3.49)  is  derived  by  the  method 
similar  to  that  in  (3.47).  The  inequality  corresponding  to 
the  second  inequality  in  (3.47)  used  to  derive  (3.49a)  holds 
with  a  large  margin  at  P.  Furthermore,  since  Y2>Y^+Y3»  the 
inequality  corresponding  to  the  second  inequality  in  (3.47) 
used  to  derive  (3.49b)  also  holds  with  a  large  margin. 

Later,  we  shall  evaluate  the  size  of  attractor  by  direct 
integration  of  equation  (3).  On  the  other  hand,  the  suffi¬ 
cient  condition  (3.45)  for  the  convergence  of  £(t)  to  9^  is 
obviously  too  strong.  The  maximum  value  of  |vAQ|  satisfying 
(3.45)  is 
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!  vA0  !  max  =  '2W<2W/U’ 
when  p3  =  |  vAQ  | /(2y2+Y3  )  and  p4  =  |  vA^  |  /  (  2y-^+Y 3  )  .  Since 
s(|vAOlmax)  <  0  ((3‘7))’  lvAolmax<  lvAol-  Hence’  if  A()2>V» 

A 

the  set  T  is  empty. 

Equation  (6)  has  a  solution  which  goes  to  infinity  as 
t  -*■ 00  on  the  x^-axis.  Hence,  there  is  no  natural  boundary  of 
the  trajectories  of  (6)  as  in  the  case  of  equation  (3).  On 
the  other  hand,  the  numerical  results  in  the  next  chapter 
suggest  that  the  pseudo-chaotic  solutions  of  (6)  are  bounded 
for  certain  values  of  parameters.  If  such  trajectories  are 

A 

contained  in  T  ((3.46)),  there  exists  a  set  0.  For  other 
values  of  parameters,  the  limit  set  of  the  pseudo-chaotic 
solutions  appears  to  be  unbounded.  If  the  limit  set  is 
unbounded,  it  is  not  contained  in  T  in  the  Y-space.  Then, 
the  sufficient  condition  in  (S.2.3)  for  the  existence  of  0 
is  too  strong,  or  there  does  not  exist  a  set  0. 

3.3.4  Local  Properties  of  Original  Systems 

In  the  previous  sections,  we  did  not  prove  the 
existence  of  0  for  system  (3')  or  (6').  In  this  section,  we 
study  the  local  properties  of  these  systems  about  certain 
sets  and  the  stablity  of  Eg ,  which  will  provide  some  insight 
on  the  convergence  of  the  trajectories  to  E. 

If  e  =  |  vAg  |  ^-y-^Y2  >  0  {3'),  then  the  equilibrium  set 

consists  of  Pqy =  {9^}  and 

EQ3={Y:  (r1}r2,r3)  =  ({(Y3/Y1)els(|vAj-ei5)}^  IvAgl-e5)}*5,  e5). 
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(3.55) 


sin(01+02)  =  sin(0.j-02-03)  =  1}, 

where  vAg  <  0  is  assumed.  The  characteristic  equation  of 
dF/dY  for  (3')  at  any  point  Py  EQ3  is  p(X)  = p  (X)p(X)  =  0 , 
where 

p(x)  =  x{x2+(y1+y2+y3)x+y3(y1+y2MvA0I 

/( |vA0|-eJ5)-2(Y3/Y1)|vAol®?5^>  (3.56) 

o  yv  9 

and  is  defined  in  (3.3b).  If  Y2  > Yi+Y3  and  Ag  >  Ag  , 
equation  (3.56)  has  a  zero  root  and  two  stable  roots.  As 

mentioned  earlier,  equation  (3.3b)  has  a  pair  of  complex 

.  •  2  '•l 

roots  with  positive  real  parts  if  y2  >  Yj+Yg  and  A0  >  Ag  • 

Hence,  all  the  eigenvectors  corresponding  to  the  unstable 

roots  of  (3.5b)  lie  on  Eg  at  Py  where  ©3 =  0.  Thus,  there  is 

no  unstable  manifold  which  is  transverse  to  E„  at  its 

corresponding  Py. 

For  system  (3'),  the  intersection  of  all  the  E  '  s  is 

0 

A 

the  origin  PQy.  The  characteristic  equation  of  dF/dY  at  PQy 
is  Pg(X)=Pg(X)  =0,  where  Pg  is  defined  in  (3.3a).  We  know 
that  (3.3a)  has  one  positive  and  two  negative  real  roots  if 
|  vAg  |  >  y-j_Y2  •  Hence,  there  is  a  two-dimensional  space  ESy+ 

spanned  by  the  eigenvectors  corresponding  to  the  unstable 
eigenvalues  at  Pgy-  As  mentioned  in  Section  3.1.1,  the 
eigenspace  in  EQ  corresponding  to  the  unstable  eigenvalues 
of  dF/dX  at  0^  is  one  dimensional.  Therefore,  ESy+  is 
transverse  to  all  the  Fg's  at  Pgy* 

We  shall  show  that  ESy+  lies  in  E  at  Pgy»  which 
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implies  that  the  trajectories  tend  to  approach  E  about  Pgy* 
The  positive  root  of  (3.3a)  is 

XY+  =  [{(y1-y2)2  +  4|vA0  |  2}1s-(y1+Y2)]/2, 

and 

ESy+  =  {Y  :  Y=  (-vAQn1,  -vAQn2,  ^Y+h2,  A^r^,  0,  0)T  ;  R}. 

On  the  other  hand,  in  E, 

yly2~ZlZ2  =  °»  (yly2+ZlZ2)y3+(zly2-ylZ2)z3  =  (3.57) 

If  Y  e  ESy+,  Y  satisfies  (3.57)  for  any  (n^ru).  Hence, 

ESy+  c  I . 

The  matrix  dF/dY  at  PQY  has  a  block  diagonal 
representation.  The  submatrix  corresponding  to  the  (y3,z3)- 
space  is  diagonal  and  has  a  repeated  eigenvalue  -Y3-  This 
implies  that  03(t)  tends  to  a  constant  about  PgY»  On  the 
other  hand,  ©.^(t)  and  02(t)  tend  to  vary  about  PqY>  since 
ESy+  does  not  coincide  with  the  (y^,z^)  or  (y^ , )-space . 

In  the  above  sense.  Eg  is  unstable  about  PqY» 

Finally,  we  examine  the  extent  of  the  unstable 
manifolds  tangent  to  ESY+  at  the  origin,  along  the  (y3,z3)- 
direction.  The  characteristic  equation  of  the  linearized 
vector  field  about  Eg  (  =  {Y  :  r^  =  r^  =  0 } )  is 

p( A )  ^  { A2+(y1+Y2 5 A+(y1Y2~ I vAq | 2 )+r32 }2  =  0.  (3.58) 

2 

Suppose  that  |  vAg  |  >  y-j_Y2  »  then  (3.58)  has  unstable  roots 

2  2 

for  r3  <  |  vAg  |  -y^Y2>  ajid  ^as  only  stable  roots  for 
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r32  >  |vA0|2-Yiy2. 


If  | vAg |  > in  (6'),  the  equilibrium  set  consists  of 


P0Y  and 


EQg  =  {Y  :  (r^r^r^  =  ((Y2Y3)^,  { ( | vAQ | -Y^Y^  >  {( | vAQ  > ; 


sin(29^)  =  sin^-O^Og)  -  1}, 


(3.59) 


where  vAQ  <  0  is  assumed.  The  characteristic  equation  of 
dF/dY  for  (6')  at  any  point  Py  in  EQg  is  p( X )  = p( A )p ( X )  =  0 , 
where 


p(A) ^  A{A2+A(y1+y2+y3+!vA0[ )+2(y2+y3) |vAq|}  =  0, 


(3.60) 


and  p  is  defined  in  (3.32b).  Equation  (3.60)  does  not  have 
unstable  roots.  Hence,  the  eigenvectors  corresponding  to 
the  unstable  roots  of  (3.32b)  lie  in  EQ  at  Py  where  0g =  0. 
Thus,  there  is  no  unstable  manifold  which  is  transverse  to 
at  its  corresponding  Py.  At  the  origin,  the  characteris- 

A 

tic  equation  of  dF/dY  has  only  one  unstable  root.  The 
eigenspace  corresponding  to  this  root  coincides  with  the 
one-dimensional  space  =  {Y  :  r2  =  r^  =  0,  cos(20^)  =  0, 
sin(20^)  =  1},  which  lies  in  Eg  for  all  9. 

For  system  (6'),  the  intersection  of  all  the  Eg's  is 
the  surface  S2  =  {Y  :  r2  =  r3  =  0 } .  The  characteristic  equation 
of  the  linearized  vector  field  about  S2  is 


p(  A )  ^  {A2  +  (Y2+Y3)*  +  Y2Y3"rl2^2  =  °* 


(3.61) 


The  above  equation  has  a  repeated  unstable  root: 


*  / 


Ay+  =  [{(Y2-Y3)2  +  4r12}is-(Y2+Y3)]/2, 

2  2 

if  >Y2Y3*  Hence,  if  >Y2Y3»  S2  is  unstable. 

Moreover,  02(t)  and  0g(t)  tend  to  vary  about  S2,  since  the 
submatrix  of  the  linearized  vector  field  corresponding  to 
the  (y2 , z2 ,y3 , )-space  is  not  diagonal  except  for  the 
origin.  In  this  sense,  Eg  is  unstable  about  S2« 

We  have  proved  that  for  (3'),  EQ  is  unstable  about  the 
origin  and  for  (6'),  Zg  is  unstable  about  Sj.  In  (6'),  the 
trajectories  always  approach  the  origin  before  approaching 
S2«  Hence,  the  trajectories  of  (3*)  and  (6')  may  not 
converge  to  one  of  Eg’s  if  the  trajectories  get  arbitrarily 
close  to  the  origin  without  converging  to  the  origin  as 
t  -*■  °°,  that  is  , 


lim  sup  |Y(t)|  >  0, 
t-***  T<t 

lim  inf  |Y(t)||  =  0, 

t-H*>  T<t 


(3.62) 


where  Y(t)  =  y(t,Yg).  The  question  is  whether  the  nonconver¬ 
gence  of  the  trajectories  to  one  of  Eg's  implies  the  noncon¬ 
vergence  of  the  trajectories  to  I.  Intuitively,  it  is  less 
likely  that  there  exists  a  trajectory  which  converges  to  the 
union  of  more  than  one  Zg's  as  t -*■«>,  but  it  is  difficult  to 
prove  it. 

3.3.5  Energy  of  Original  Systems 

Let  I =  { 2 | | 2+ j A2 | 2+ | Ag | 2 } /2  in  system  (3').  We  call 


I  an  energy  function  of  (3')  for  convenience.  We  shall  see 
the  difference  between  the  behaviors  of  I  for  the  original 
and  reduced  systems.  Taking  the  time  derivative  of  I, 

i  =  -[2Y1{  |A1|-(3v/4y1)  |vAq1  JsinO^+Qj)  |  |a2  |  }2 

+{y2-(3v|vA0| |sin(ei+e2) j)2/(8Y1>}|A2|2+Y3|A3|2],  (3.63) 

where  v =  sgnfvA^sinC 0^+0 2 ) ) •  Equation  (3.63)  is  sketched  in 

2 

Figure  3.7,  where  E  is  a  constant  I  surface.  Since  |vAq| 

>  y1y2’  t^le  seconc^  term  of  the  right  hand  side  of  (3.62)  can 


I  A*  I  |A3 1 


Figure  3.7  Dynamics  of  Energy  of  (3') 

be  negative.  Hence,  we  have  an  elliptic  cone  C  :  I  =  0,  in 
which  the  energy  I  increases,  and  outside  which  I  decreases. 
The  line  of  foci  and  the  radii  of  the  cone  vary  as  9]_+92 
varies.  We  assume  that  vAq  <  0.  Then,  if 


sin(01  +  02)  <  (8^/3)  (y1Y2  i55/  |  vAQ  |  , 

cone  C  does  not  exist  and  the  energy  I  decreases.  For  (3) 
or  on  Z,  sin(0^+©2)  =  ±1  and  the  cone  C  is  fixed. 

Taking  the  time  derivative  of  I  for  (6'), 

I =  -{2(Y1-v|vAo||sin201|)|A1|2+Y2|A2|2+Y3|A3|2},  (3.64) 

where  v = -sgn(vAQsin201) .  Since  |vAQ|  >Y^  for  chaotic 
motion,  the  first  term  of  the  right  hand  side  of  (3.64)  can 
be  negative.  Hence,  we  also  have  an  elliptic  cone  C  as  in 
the  case  of  (3M.  The  line  of  foci  is  fixed  at  the  |A^|- 
axis,  and  the  radii  of  cone  vary  as  0^  varies.  If  sin20^ 
<Y1/|vAq|,  I  always  decreases.  For  (6')  or  on  £,  sin201  =  ±1 
and  the  cone  C  fixed. 
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CHAPTER  4 


NUMERICAL  EXPERIMENTS 


In  this  chapter,  we  present  some  numerical  results  for 
equations  (3'), (3), (6')  and  (6).  An  attempt  will  be  made  to 
correlate  some  of  the  numerical  results  with  the  analytical 
results  given  in  the  previous  chapters.  Equation  (3)  has 
already  been  studied  numerically  by  Pikovskii  et  al.  [12], 

In  the  calculations,  we  use  the  normalized  time  yt,  where 
y  =  max{y1,Y2 »Y3)  for  (3')  and  (3),  and  y  =  max{y2,y3}  for 
(6')  and  (6).  In  what  follows,  we  denote  the  normalized 
time  by  t  also.  By  the  substitutions: 


A^ / y  **■  A^ ,  i  =  0,1, 2, 3, 

x./y  xi *  Yj/Y  *  Yj.*  i  =  1,2,3,  , 


(4.1) 


we  obtain  a  system  described  by  the  same  equation  except  for 
Y  =  1 .  Any  equations  derived  in  the  previous  chapters  are 
not  changed  by  this  substitution. 


4.1  Behavior  of  Original  Systems  (3')  and  (6') 

We  fix  the  values  of  parameters  as  follows  :  y^  =  Y3  = 
0.25,  y2 =  1.  These  are  the  values  used  by  Pikovskii  et  al. 
[12]  for  the  model  (1.4)  describing  the  parametric  interac¬ 
tion  of  a  whistler  with  ion  sound  and  plasma  oscillations 
near  the  lower  hybrid  resonance  combined  with  the  three-wave 
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interaction  involving  another  plasma  wave  synchronous  to  the 
parametrically  excited  pair.  For  these  values,  the  condi¬ 
tion  Y2  >  is  satisfied.  Hence,  the  stability  of  the 

equilibrium  points  is  as  follows  :  For  | vAQ |  <  |vAQ|  =0.5, 

PQ  is  stable  and  does  not  exist  ;  for  0.5  <  |vAg|  <  j vAQ | 

=  1.259,  PQ  is  unstable  and  P^  is  stable  ;  and  for  |vAQ| 

>  1.259,  Pq  and  P^  are  unstable. 

The  condition  (3.45)  is  violated  for  | vAQ  |  >  1.259  . 

On  the  other  hand,  numerical  results  show  that  the  phases 
are  nearly  locked  for  j vAg |  <3.0  at  least  for  150  time  units 
The  plots  of  cos(01(t)-02(t)-02(t) )  and  (Re(A^( t ) ) /  |  vAQ  |  , 
Im(A^(t ) ) / | vAq I )  are  shown  in  Figures  4.1  and  4.2,  respect¬ 
ively.  As  | vAg |  increases  from  3.0,  the  nearly  phase-locked 
state  appears  to  be  unstable.  At  |vAQ|  =3.2,  an  abrupt 
phase  change  occurs  occasionally  in  an  unpredictable  manner 
(Fig. 4 . 1(a) ,4 . 2 (a) , (a* ) ) .  In  Figure  4.2(a'),  we  do  not 
observe  any  abrupt  phase  change  for  150  time  units,  while  in 
(a)  (corresponding  to  Fig. 4. 1(a)),  it  is  observed  twice. 

From  this  observation,  it  seems  to  be  difficult  to  determine 
numerically  the  critical  value  of  | vAQ |  at  which  such  a 
phase  change  takes  place.  For,  even  if  no  abrupt  phase 
change  is  observed  over  a  long  time  interval,  we  do  not  know 
whether  it  will  occur  at  the  next  moment.  At  |vAq|  =3.5, 
the  abrupt  phase  change  occurs  more  often  (Fig. 4. 1(b), 
4.2(b)).  At  | vAq I  =4.0,  the  phase  varies  in  a  complicated 
manner  most  of  the  time  (Fig. 4 . 1(c) ,4 . 2(c) ) .  At  |vAq|  =  8.0, 
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the  phase  varies  periodically  with  time  ( Fig. 4 . 1 (d)  ,4 . 2 (d ) ) . 
Thus,  for  | vAq |  >3.5,  the  reduced  system  (3)  does  not 
represent  the  original  system  (3')  most  of  the  time.  For 
|  vAq  |  <3.2,  system  (3)  may  represent  (3')  most  of  the  time, 
but  not  asymptotically. 

It  is  observed  that  the  abrupt  phase  changes  occur 

only  about  Eg  near  the  origin.  This  is  consistent  with  the 

result  in  Section  3.3.4  where  it  was  shown  that  there  exists 

2  2 

an  unstable  manifold  at  Eg  if  r3  <  |  vAQ  |  “Y^Y2’  The  extent 
of  this  manifold  is  increased  as  | vAg |  increases,  which  may 
be  one  of  the  reasons  why  the  abrupt  phase  changes  occur 
more  often  as  |vAQ|  increases.  According  to  the  observation, 
the  hypothesis  (3.62)  appears  to  be  satisfied,  that  is,  the 
trajectories  always  return  to  an  arbitrary  neighborhood  of 
the  origin.  Furthermore,  the  numerical  results  show  that 

2  /s  2 

as  Ag  decreases  to  AQ  ((3.8)),  the  trajectories  are 
trapped  in  a  smaller  neighborhood  of  E,  but  they  do  not 
converge  to  E  as  t-*-00. 

For  system  (6'),  we  fix  the  values  of  parameters  as 
follows  :  =  y2  =  0.4  and  Y3  =  1*  These  values  are  reasonable 

for  the  model  such  that  plasma  wave  1  generated  by  an 
external  wave  decays  into  another  plasma  wave  2  and  an  ion 
acoustic  wave  3.  The  numerical  results  show  that  for  a  wide 
range  of  parameter  values,  the  phase  becomes  locked  rapidly 
and  remain  locked  (within  the  accuracy  of  computation) . 

The  phase  is  locked  even  when  condition  (3.47)  is  not 
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As  mentioned  in  Section  3.3.4,  the  trajectories 
approach  the  origin  before  approaching  From  the  obser¬ 

vation,  equation  (3.62)  appears  to  be  satisfied  for  certain 
values  of  | vAg | .  Therefore,  the  trajectories  may  return  to 
an  arbitrary  neighborhood  of  S2  repeatedly.  According  to 
the  analysis  in  Section  3.3.4,  EQ  is  unstable  about  S„. 
Hence,  the  abrupt  phase  changes  may  occur  about  S2>  But, 
actually,  we  did  not  observe  any  abrupt  phase  changes  for 
( 3 ' )  for  a  wide  range  of  parameter  values.  This  result  may 
be  explained  as  follows.  Since  the  trajectories  approach  E 
exponentially  about  the  origin  (Section  2.3),  therefore  as 
a  trajectory  gets  closer  to  S2 ,  it  is  closer  to  E.  From  the 
continuity  of  the  vector  field,  the  vector  F(Y)  tends  to 
parallel  to  E0  as  Y  approaches  E  except  in  the  neighborhood 
of  a  certain  set  containing  S2  and  EQg.  The  volume  of  such 
a  neighborhood  is  smaller,  as  a  smaller  neighborhood  of  E 
is  considered.  Hence,  it  is  less  likely  that  the  trajectory 
moves  transverse  to  Eg  as  it  gets  closer  to  E.  That  is,  it 
is  less  likely  that  about  S2,  the  trajectory  moves  across 
Eq's  and  the  abrupt  phase  change  occurs. 

4.2  Behavior  of  Reduced  Systems  (3)  and  (6) 

In  what  follows,  we  consider  only  the  reduced  system 
( 3 )  and  ( 6 ) .  For  ( 3 ) ,  we  assume  that  1 . 2  59  <  | vAg |  <3.2  for 
which  the  system  may  have  a  chaotic  solution,  and  describes 


(3')  approximately  most  of  the  time,  but  not  asymptotically 
as  t -*•<*>.  For  (6),  we  assume  that  |  vAq  |  >0.4  for  which  (6) 
may  have  a  chaotic  solution  and  describes  the  asymptotic 
behavior  of  (S’). 

4.2.1  System  (3) 

The  projection  of  the  trajectory  onto  the  (x^,Xg)- 
plane  is  shown  in  Figure  4.3,  At  each  value  of  vA^ ,  there 
appears  a  pseudo-chaotic  attractor.  According  to  Pikovskii 
et  al .  [12],  periodic  solutions  appear  for  | vAq |  £3.35. 

But,  as  mentioned  earlier,  equation  (3)  does  not  describe 
equation  (3')  approximately  most  of  the  time  for  such  values 
of  | vAQ ! . 

The  attractors  in  the  Figure  4 . 3  resemble  two-dimen¬ 
sional  surfaces.  But,  their  cross  sections  have  a  complica¬ 
ted  structure.  Let  us  consider  a  connected  set  Q  on  a  plane 
PL  which  is  transverse  to  the  trajectories  as  shown  in 
Figure  4.4.  In  equation  (3),  the  phase  volume  shrinks 
uniformly,  since 

3x..  3x„  3x 

3^  *5^  *3^  -  -<WV  < 

Hence,  the  area  of  Q  shrinks  exponentially  to  zero  as  t-*00. 
Simultaneously,  as  shown  in  Figure  4.3,  Q  is  stretched  in 
one  direction  while  moving  around  the  nontrivial  equilibrium 
point  Py  ( Fig . 4 . 3 (a) ) .  Hence,  compression  of  Q  must  take 
place  along  another  direction.  In  Figure  4.3(a),  Q  is  split 
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into  two  parts  at  the  origin  and  then  folded.  In  (b),  Q  is 
also  bent  along  those  portions  of  the  trajectories  which  are 
close  to  Py .  If  we  consider  a  first  return  mapping  from  PL 
into  PL,  Q  is  finally  mapped  into  a  set  exhibiting  a  Cantor- 
set  like  structure  after  stretching,  splitting  and  bending 
infinitely  many  times. 


4.2.2  System  (6) 

Here,  a  detailed  numerical  study  for  system  (6)  is 
made.  The  projection  of  the  trajectories  onto  the  (x^,x3)- 
plane  is  shown  in  Figure  4.5.  In  equation  (6),  the  phase 
volume  also  shrinks  uniformly  if  |vAQ|  < 3  =  1 . 8 ,  since 


3x. 

3x~ 

3x_ 

_ 1 

+ _ 2 

+  — * 

3x^ 

3x2 

3x3 

f vAq 1  - ( y 1+y 2  +y 3 )  <  0. 


(4.2) 
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Figure  4.5  Projections 


't 


As  =  |  v  A  gl~ =  0.2,  the  solution  is  pseudo-chaotic 
(see  Fig. 4. 5(a)).  A  set  Q  as  defined  in  the  previous 
section  is  split  at  the  origin  and  is  bent  along  the  closest 
trajectory  to  PV  .  The  closure  of  the  attractor  seems  to 
contain  the  x^-axis,  which  means  that  the  attractor  is 
unbounded.  As  |  vAq  |  increases,  the  trajectory  tends  to 
leave  P^  more  rapidly  in  a  spiral  manner.  Actually,  the 
ratio  r  of  the  frequency  and  the  growth  rate  of  the  spiral 
trajectory  at  decreases  as  |vAq|  increases,  since  from 

equation  (3.33b), 

dr2  (r2-3)2(r2+9)(Y1+Y2+Y3“lvAol )2^2( lvAol“Yl)+Y2+T3^ 

d|vAQ|  4(l+r2)( | vAq |-y1)2y2y3 

Thus,  at  y^  =  0 • 4 ,  the  pseudo-chaotic  attractor  is  far  from 
the  origin  (Fig4.5(b)).  Hence,  the  set  Q  is  only  bent  along 
the  trajectory  closest  to  .  At  y^  -  0  •  488  'v  0 . 4896  (Fig. 
4.5(d)),  two  pseudo-chaotic  attractors  appear.  They  are 
linked,  but  disconnected  with  each  other.  In  the  figure, 
one  of  them  is  omitted.  The  omitted  one  is  symmetric  to  the 
given  one  with  respect  to  a  180°  rotation  around  the  x2~axis. 
At  Y^  =  0.4875,  these  two  pseudo-chaotic  attractors  have  a 
small  intersection  so  that  the  trajectory  moves  from  one  to 
the  other  infrequently  as  shown  in  Figure  4.5(c).  At  y^ = 
0.4896  ^  0.4899  ,  a  three-loop  pseudo-chaotic  solution  appears 
(  Fig  .  4 . 5  (e ) )  .  At  Y]_  =  0.49,  we  have  a  three-loop  periodic 
solution  ( Fig . 4 . 5 (f ) ) ,  where  the  trajectory  in  a  transient 
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state  should  be  neglected.  At  =  0 . 491  %  0 . 495  ,  a  pseudo- 
chaotic  attractor  of  the  type  in  Figure  4.5(d)  appears  again 
( Fig . 4 . 5 ( g) ) .  At  0.4955,  an  eight-loop  solution  appears 

(Fig. 4 . 5 (h) ) .  At  -  0 . 496  v  0 . 497  and  0.498  ,  we  have  a 
four-loop  pseudo-chaotic  attractor  (Fig.4.5(i))  and  a  stable 
four-loop  periodic  solution  ( Fig  .  4 . 5 ( j ) )  ,  respectively. 
Furthermore,  at  7^  =  0.49  8^0.5  and  0.51,  we  have  a  two-loop 
pseudo-chaotic  attractor  (Fig.4.5(k))  and  a  two-loop 
periodic  solution  ( Fig . 4 . 5 ( 1 ) ) ,  respectively.  At  y^-0.52 
v  0.56,  there  is  a  single-loop  periodic  solution  (Fig. 4. 5 
(m)).  This  solution  is  not  accompanied  by  a  single-loop 
pseudo-chaotic  solution.  As  y|  increases,  this  asymmetrical 
single-loop  periodic  solution  gradually  shifts,  and  finally, 
a  symmetrical  one  appears  at  y^>0.6. 

4.2.3  Sizes  of  Attractors 

In  Figure  3.6,  P=  (0,  -S^vAg,  -2vAQ).  Hence,  for  the 
values  of  vAQ  given  in  Figures  4. 3(a), (b)  and  (c),  P= 

( 0 , 2  .  53 , 3 . 56  ) , (  0 , 7 . 07 , 5 . 0  )  and  (0,9.05,6.4),  respectively. 
From  the  numerical  experiments,  the  maximum  values  of  | x  2  ^  t ) | 
and  x^(t)  with  respect  to  t  are  approximately  (2. 0,2.0), 

(4. 0,4.0)  and  (4. 5, 5.0),  respectively.  Hence,  the  sizes  of 

A 

attractors  are  smaller  than  A  as  predicted  in  Section  3.3.3. 

Rewriting  equation  (3.46)  for  y2  =  0.4,  y3  =  1, 

| A2  j  2  + I A3 | 2  <  1.4y1{l-|vA0|2/(y1  +  1.4)2}.  (4.4) 

Hence,  for  fixed  Y-[  =  ( vA  q  I  —  Y  »  the  size  of  f  depends  on  Yj_  • 


Assume  that  y^ =  0.489,  for  instance.  Then,  from  Figure  4.5 
(d),  the  maximum  value  of  | A2  I  is  approximately  5.49. 

The  right  hard  side  of  (4.4)  has  a  maximum  value  of  2.546 
in  the  limit  y^-*-°°.  Hence,  equation  (4.4)  is  not  satisfied. 
This  result  implies  that  the  sufficient  condition  (3.46)  for 
the  convergence  of  £(t)  to  9^  is  not  applicable. 

4.2.4  Periodic  Solutions  and  Bifurcations 

The  asy’  itric  periodic  solution  in  Figure  4.5(m) 
successively  bifurcates  as  y^  decreases.  The  observed 
transition  for  y^ e  [0.489,0.5)  resembles  that  of  the  well 
known  one-dimensional  mapping  [21,22] 

x(k+l)  =  G(x(k))  I  rx(k) ( l-x(k ) )  (4.5) 

which  has  properties  shared  by  more  general  models  [28]. 
Actually,  we  obtain  a  first  return  mapping  of  the  trajectory 
for  y^  =  0.489  in  Figure  5.1(d),  which  seems  to  satisfy  the 
conditions  for  the  models  in  reference  [28]. 

Figure  4.6  is  a  schematic  diagram  of  bifurcation 
phenomena  of  system  (6).  Here,  the  circles  represent 
actually  observed  periodic  solutions,  but  the  branching  is 
conjectured  from  that  of  (4.5).  The  broken  lines  are 
unstable  closed  orbits  and  the  solid  lines  are  stable  closed 
orbits.  From  Figure  4.6,  the  bifurcation  phenomena  of  (6) 
for  y^  e  [0.489,0.55)  seem  to  be  explainable  by  a  one-dimen¬ 
sional  mapping  of  type  (4.5),  although  it  is  a  rough 
approximation . 
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Figure  4.6  Schematic  Diagram  of  Bifurcation  of  (6) 

According  to  [21,22,28],  the  "chaos"  appears  when  the 
bifurcation  parameter  exceeds  an  accumulation  point  of  the 
sequences  of  bifurcation  values,  which  correspond  to  Figures 
4.5(d)  and  (g).  For  other  values  of  parameters,  stable 
periodic  orbits  appear.  Hence,  the  three-loop  pseudo- 
chaotic  solutions  for  -  0 . 4896  'v  0 . 4899  seem  to  be  3*2n 
time-loop  stable  periodic  orbits,  and  the  two  and  four-loop 
pseudo-chaotic  solutions  in  Figures  4.5(i)  and  (k)  seem  to 
be  converging  to  a  stable  periodic  orbit. 

4.3  Simple  Models  Describing  Attractors 

In  this  section,  we  consider  the  existence  and  the 
transition  of  the  attractors  shown  in  Figures  4.3  and  4.5. 

We  consider  only  the  attractors  in  which  there  occurs  the 
splitting  and  bending  mentioned  earlier,  since  such  mechan- 


100 


isms  may  lead  to  chaotic  behavior.  It  is  too  difficult  to 
verify  mathematically  the  existence  and  transition  of  the 
attractors.  Therefore,  we  introduce  a  simple  two-dimen¬ 
sional  model  whose  behavior  is  similar  to  that  of  (3)  and 
(6).  This  model  is  only  a  schematic  one  ,  but  is  useful  to 
explain  the  behavior  of  the  attractors  and  to  understand  the 
relation  between  the  trajectories  in  Figures  4.3  and  4.5 
and  the  first  return  mappings  shown  in  Chapter  5. 

We  consider  half  planes  and  K2  containing  the 
unstable  foci  F^  and  F^,  respectively  (Fig. 4. 7  and  4.8). 

The  trajectories  move  away  from  F^  (resp.,  F 2)  in  a  spiral 
manner  on  (resp.,  K2)  and  jump  to  the  line  L2  (resp.,  L^) 
parallel  to  the  x2~axis  in  Figure  4.7  and  to  the  x3~axis  in 
Figure  4.8.  The  basic  assumptions  are  as  follows: 

(A. 4.1)  Foci  F^  and  F2  move  upward  slowly  as  | vAQ | 
increases  so  that  it  crosses  and  L2,  respectively; 
(A. 4. 2)  The  ratio  r  of  the  frequency  and  the  growth 
rate  of  the  spiral  trajectory  about  the  nontrivial 
equilibrium  points  decreases  rapidly  as  |vA0| 
increases ; 

(A. 4. 3)  Suppose  that  t^  <  t2 ,  X^  =  X(t^)  e  ,  i  =  1,2 
and  X(t)  t  L^,L2  for  all  t«  (t^,t2).  Then,  |  X2||  /  |  X^|| 
decreases  monotonically  as  fX^fl  increases,  and  if 
|X2||/|X^||  >1  for  sufficiently  small  1  X^fl  ,  fx^/lx^  =  1 
for  a  certain  value  of  ||X^||.  Here,  suffices  1  and  2 
can  be  exchanged. 
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Actually  for  (6),  the  X2  and  x,  components  of  increase 

perpendicularly  to  the  x^-axis  as  |vAg|  increases.  For  (3), 
behaves  in  a  similar  way  for  sufficiently  large  |vAg|. 
Numerical  results  and  equation  (4.3)  show  that  (A. 4. 2)  is 
reasonable  for  the  models  describing  (3)  and  (6). 

Assumption  (A. 4. 3)  implies  that  the  trajectory  does  not 
diverge  as  it  bounces  between  and  Lj  without  passing 
through  F^O  or  F2O.  We  shall  introduce  other  assumptions 
for  each  of  (3)  and  (6)  later. 

This  model  resembles  a  "universal  circuit"  used  in 
reference  [29],  but  is  not  obtained  by  taking  limits  or 
choosing  specific  values  of  parameters  in  (3)  or  (6).  That 
is,  introducing  such  a  model  is  not  justified  mathematically 
or  from  the  physical  point  of  view.  We  use  this  model  to 
avoid  the  difficulties  in  studying  equations  (3)  and  (6) 
directly. 

For  the  model  describing  (3),  we  assume  that  the 
trajectory  converging  to  the  origin  as  t -*■ encircles  F^ 
(resp.  ,  F2>  and  hits  Lj  (resp.,  L^)  at  B2  (resp.,  B1)  near 
F2  (resp.,  F-^).  Then,  by  decreasing  r  according  to  assump¬ 
tion  (A.4.2),  we  obtain  the  transition  as  shown  in  Figure 
4.7.  Here,  the  broken  and  dotted  lines  are  trajectories, 
and  at  (resp.,  Gj),  they  are  tangent  to  (resp.,  Lj). 
The  shadowed  regions  in  Figure  4.7  are  obviously  local 
attractors.  In  (a)  and  (b) ,  they  are  also  global  attractors 
by  assumption  (A. 4. 3).  Figures  4.7(a)-(c)  correspond  to 
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Figures  4.3(a)-(c),  respectively.  In  Figure  4.7(a),  the 


line  segment  (corresponding  to  the  set  Q  in  Fig. 4. 4)  is 

split  at  the  origin.  In  (b),  a-*-so  bent  at  G^  and  G^. 

By  assumption  (A. 4.1)  and  (A. 4. 2),  as  |vAg|  increases, 
points  A^  and  in  (a)  move  upward  more  rapidly  than  G^, 
and  moves  downward  first,  hits  the  origin  and  then  moves 
upward.  We  assume  that  hits  the  origin  after  B^  and  G^ 
coincide  at  some  value  of  | vA_ | .  When  B.  coincides  with  G., 

0  1  i 

the  attractor  of  type  (a)  disappears.  When  B^  is  above  G^, 
there  appears  as  an  attractor  of  type  (b)  or  a  multi-loop 
attractor  of  type  (c).  After  hits  the  origin,  the 
situation  is  the  same  as  that  of  the  model  describing  (6). 

The  model  describing  (6)  is  shown  in  Figure  4.8.  We 
assume  that  the  x^-axis  is  a  trajectory  and  x^t)  = 
x1(0)exp(y^t) ,  y^>0.  Hence,  in  (a),  the  attractor  is 
unbounded  and  there  exist  various  types  of  trajectory  behav¬ 
ior  as  shown  in  the  figure.  The  line  segment  is  split 

at  the  origin  and  bent  at  G^  and  G2.  This  case  corresponds 
to  Figure  4.5(a). 

As  | vA g ]  increases,  F^  moves  upward  and  moves 
downward,  hits  the  origin  and  then  moves  upward.  Assume 
that  all  the  points  are  located  as  shown  in  (b).  Such  an 
attractor  corresponds  to  that  in  Figure  4.5(b).  The 
attractor  is  bounded,  and  the  line  segment  P^Q^  is  bent  at 

G.  . 

1 

As  |vAq|  increases  further,  A^,  B^  and  approach  G^. 


*  / 
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Assume  that  and  approach  G^  more  rapidly  than  B^. 

Then,  A^  and  exchange  their  positions  at  a  certain  value 
of  | vAq | ,  before  A^  and  hit  G^.  Then,  we  have  an 
attractor  of  type  (c),  where  another  attractor  obtained  by 
180°  rotation  about  the  X2~axis  is  omitted.  Hence,  the  line 
segment  E-^H^  contains  B-^C^  (which  is  omitted  in  the  figure), 
i.e.,  the  attractors  are  interlinked  and  disconnected.  When 
A^  is  slightly  above  B^,  we  have  an  attractor  of  type  (b) 
corresponding  to  Figure  4.5(c). 

When  A^  is  slightly  above  G^,  we  have  a  two-loop 
attractor  as  shown  in  (d).  Moreover,  in  the  attractor  of 
type  (c),  we  may  have  multi-loop  attractors.  Assume  that 
the  attractor  is  n  -loop  and  encircles  G„  n  times  on  K„ . 
Obviously,  for  the  existence  of  an  attractor  in  which  there 
occurs  a  bending,  A^  is  above  B£  if  np-ne  is  odd,  and  A£  is 
below  if  np-ng  is  even  (Figure  4.9).  Practically,  it  is 


Figure  4.9  Multi-loop  Attractors 


likely  that  there  exists  a  stable  periodic  solution  in  such 
a  small-width  multi-loop  attractor.  In  the  numerical 
experiment,  we  did  not  find  any  evidence  for  the  existence 
of  such  a  multi-loop  attractor. 

If  C2  is  above  G 2  for  certain  value  of  |vAg|  as  shown 
in  (e),  the  attractor  of  type  (c)  does  not  exist.  Finally, 
assume  that  is  above  as  shown  in  (f).  From  assumption 
(A. 4. 3),  there  exists  a  point  where  ! X2 | / fl X, |  =  1.  Then, 
the  mapping  of  C2D2  (resp.,  C^D^)  into  C-^D.^  (resp.,  C2D2)  is 
a  contraction  mapping  and  we  have  a  single-loop  symmetric 
stable  periodic  solution  corresponding  to  Figure  4.5(1). 
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CHAPTER  5 


FIRST  RETURN  MAPPINGS  AND  STATISTICAL  PROPERTIES 


Most  of  the  analytical  works  on  the  chaotic  solutions 
of  nonlinear  systems  have  been  devoted  primarily  to  discrete¬ 
time  systems  [21,22,28,30-32].  For  continuous-time  systems 
with  dimension  greater  than  two,  an  analytical  approach  to 
the  problem  of  existence  of  chaotic  solutions  appears  to  be 
quite  formidable.  Up  to  the  present  time,  most  of  the  works 
on  such  systems  resort  to  numerical  experimentation.  The 
numerical  results  show  that  many  three-dimensional 
continuous-time  systems  have  pseudo-chaotic  solutions  lying 
in  sets  resembling  two-dimensional  surfaces  in  the  system's 
state  space.  For  such  cases,  we  obtain  a  quasi  one¬ 
dimensional  mapping  by  taking  a  first  return  mapping  (Fig. 
4.4).  In  this  chapter,  we  study  such  a  quasi  one-dimensional 
mapping  for  simplicity. 

5.1  First  Return  Mappings  of  (3)  and  (6) 

First  return  mappings  of  the  trajectories  in  Figures 
4. 3(a), (b)  and  4 . 5 ( a ) , (b ) , ( d)  are  given  in  Figures  5.1(a)- 
(e),  respectively.  For  the  trajectories  in  Figures  4.3(a) 
and  (b),  we  consider  mappings  from  T+  =  (X  :  x2  =  x2g,  x3  > x3e} 
or  T_  =  {X  :  x2  =  x~g ,  x3  >  x3g}  to  T+  or  T_ ,  where  (x*g ,x2g ,x3e ) 
are  equilibrium  points  defined  in  (3.2)  and  xtg>0,  x7g <  0 , 
i=  1,2.  Such  a  mapping  is  the  first  return  mapping  of  the 
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trajectories  X=  ( | x^( t ) | , | ( t ) | ,Xg ( t ) ) .  The  behaviors  of 

the  modified  trajectories  are  identical  to  those  of  the 

original  ones  because  of  the  system’s  symmetry  around  the 

Xg-axis.  For  Figure  5.1(c),  we  choose  T+  =  (X  : 

xx0  /(x*  -4) ,  x*  <  x.  <  4}  and  T  =  (X  :  x  =  (x, -xT  )  x„ 

2e  le  le  1  -  2  1  le  2e 

/(x~e+4),  -4<x1<x“e},  where  ^  xie  >x2e  ,x3e  ^  are  defined  i-n 
(3.31)  and  xt  >0,  xT  <0,  i =  1, 3 .  In  (d),  T={X:x_=0, 

X6  2.6  *  w 

x^  >  0 }  and  T  =  { X  :  x^  =  0 ,  x^  <  0  }  .  In  (e),  we  choose  a 
mapping  from  T+  =  {X  :  x^  =  0,  x^  >  0}  into  itself.  We  should 
note  that  T+  are  transverse  to  the  trajectories.  The 
intersection  of  the  trajectories  with  T+  corresponds  to 
those  on  P^Q^,  i =  1,2  in  Figures  4.7  and  4.8. 

The  mappings  given  in  Figures  5.1(a),  (c)  and  (e) 
resemble  one-dimensional  mappings.  In  (b)  and  (d),  we 
observe  that  the  mappings  are  more  complicated.  For  (b)  and 
(d),  there  exist  no  simple  planes  like  T+  on  which  we  have 
a  first  return  mapping  resembling  an  one-dimensional  mapping. 
This  is  due  to  the  fact  that  the  bent  Q  (Fig. 4. 4)  is  not 
completely  folded  after  Q  encircles  the  origin  from  T+  (or 
T_ )  to  T_  (or  T+ ) . 

In  Figures  5.1(a)-(c),  the  whole  interval  of  x^ 
contains  the  points  which  describe  the  trajectories 
converging  to  the  origin  as  t -*■<*>.  Near  the  origin,  the 
velocity  can  be  arbitrarily  small.  Hence,  for  these  cases, 
the  continuous-time  systems  (3)  and  (6)  cannot  be  replaced 
by  discrete-time  systems.  In  (d)  and  (e),  there  are  no  such 
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points  so  that  there  is  a  supremum  for  the  time  in  which  the 
trajectories  return  to  T+ . 

In  Figures  5 . 1  ( a '  )  -  ( e  '  )  ,  we  show  one-dimensional 
mappings  of  simple  models  in  Figures  4.7  and  4.8.  The  half 
planes  T+  and  T  correspond  to  and  ^2^2  figures  4.7 

and  4.8,  respectively.  These  mappings  are  useful  in 
conjecturing  the  behavior  of  the  first  return  mappings  of 
(3)  and  (6)  at  those  points  which  are  not  observed  in  the 
experiment . 

5.2  One-Dimensional  Mappings 

We  regard  the  first  return  mappings  in  Figures  5.1(a), 
(c)  and  (e)  as  one-dimensional  discrete-time  systems.  As 
mentioned  in  the  previous  section,  this  means  that  for  (a) 
and  (c),  the  temporal  behaviors  of  the  continuous-time 
versions  of  systems  (3)  and  (6)  are  not  considered. 

One-dimensional  discrete-time  systems  having  chaotic 
solutions  have  been  studied  extensively  in  mathematics  [28, 
30-32]  and  other  areas  [21,22].  Here,  we  present  some  of 
the  known  results  and  discuss  their  applicability  to  the 
systems  under  consideration. 

5.2.1  Period  Three  Implies  Chaos 

The  following  theorem  is  due  to  Li  and  Yorke  [31]: 

(S.5.1)  Let  J  be  an  interval  and  let  G  :  J  -*■  J  be 
continuous.  Assume  that  there  is  a  point  ae  J  for 
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2  3 

which  the  points  b  =  G(a) ,  c  =  G  (a)  and  d  =  G  (a) 

satisfy 

d  <  a  <  b  <  c  (or  d>a>b>c). 

Then, 

T1  :  for  every  k =  1,2,...,  there  is  a  periodic  point 
in  J  having  period  k. 

Furthermore , 

T2  :  there  is  an  uncountable  set  S  =  J  (contains  no 

periodic  points),  which  satisfies  the  following 
conditions : 

(A)  For  every  p,q £  S  with  p  i  q, 

(T2.1)  lim  sup  |Gn(p)  -  Gn(q) |  >  0 
n-*<>» 

and 

(T2.2)  lim  inf  | Gn(p )  -  Gn(q) |  =  0 

n-M« 

(B)  For  every  pc  S  and  periodic  point  qe  J, 

(T2.3)  lim  sup  | Gn(p)  -  Gn(q) |  =  0 
n-*°° 

If  there  is  a  periodic  point  with  period  3,  then  the 
hypothesis  of  the  theorem  are  satisfied.  Statement  (A) 
means  two  sequences  (Gn(p)}  and  (Gn(q)},  p,qe  S  can  be  made 
arbitrarily  close  to  each  other  for  sufficiently  large  n, 
but  they  do  not  converge  to  each  other  as  n -*•<*>.  Statement 
(B)  implies  that  a  sequence  (Gn(p)},  pc  S  does  not  converge 
to  any  periodic  points.  If  S  is  an  attractor,  it  is 
obviously  a  strange  attractor. 
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The  mappings  in  Figures  5 . 1 (a) ,  (c ) ,  (e  )  and  (a')-(e') 
satisfy  the  hypotheses  of  (S.5.1),  and  therefore,  there 
exists  a  S  for  each  of  them.  Theorem  (S.5.1)  does  not  give 
any  information  about  the  size  of  S .  If  there  is  a  stable 
periodic  point  Pg«  J,  then  there  is  a  neighborhood  N(pg)  of 
pg  such  that  G^(p)  e  W(pg)  implies  Gn(p)  e  N(pg)  for  all  n  >  N. 
At  the  present  time,  it  is  not  known  whether  the  following 
statement  is  true  or  false: 

(S.5.2)  The  existence  of  a  stable  periodic  point 
Pg€V,  where  V  is  an  attractor,  implies  that  the 
sequence  {Gn(p)}  converges  to  Pg  as  n  ■+•  ®  for  almost 
all  p  c  V. 

If  this  statement  is  true,  then  even  if  the  trajectories  in 
the  transient  state  appears  to  be  very  complicated,  there  is 
no  strange  attractor. 

Many  examples  of  one-dimensional  mappings  derived  from 
actual  systems  have  points  pg  at  which  dG(p)/dp=  0  as  that 
in  Figure  5.1(e).  There  may  exist  a  stable  periodic 
solution  Pg  near  Pg.  Experimental  results  show  that  there 
exists  a  complicated  solution  which  is  not  periodic  (within 
the  accuracy  of  computation).  This  does  not  necessarily 
mean  the  nonexistence  of  Pg  or  that  (S.5.2)  is  false.  For, 
even  if  there  exists  a  Pg  and  (S.5.2)  is  true,  a  small 
perturbation  in  the  system  may  shift  the  trajectory  out  of 
N(pg)  if  the  size  of  N(pg)  is  very  small.  In  other  words, 
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from  the  practical  standpoint,  the  existence  of  Pg  may  be 
ignored  for  certain  cases. 

We  consider  the  mappings  in  Figures  5. 1(b),  (d).  In 
these  cases,  the  mapping  cannot  be  approximated  by  an  one¬ 
dimensional  mapping  especially  about  Pg.  We  do  not  know  how 
such  a  structure  about  Pg  affect  the  existence  or  nonexist¬ 
ence  of  stable  periodic  orbits  in  the  three-dimensional 
space.  It  should  be  noted  that  the  mapping  in  (e)  also  has 
such  a  structure  about  Pg  which  is  not  so  explicit  as  those 
of  (b)  and  (d). 

5.2.2  Statistical  Properties  of  {Gn(p)} 

We  consider  the  following  theorem  given  in  [31]: 

(S.5.3)  Let  G  :  J -*•  J  satisfy  the  following  conditions: 

(A)  G  is  continuous. 

(B)  Except  at  one  point  p  e  J,  F  is  twice  continuouly 
differentiable . 

(C)  inf  |dG(q)/dq|  >  1. 
qe 

Then,  there  is  a  unique  function  g  :  such 

that  for  almost  all  p  £  J,  g  is  the  density  of  p. 

Here,  g  is  defined  to  be  a  density  of  p  if 

4><P»Cp,  ,p,])  =  lim  4>(p,N,[p,  ,p9]>  =  j?2g(q)dq 
N-*»  pl 

for  all  p^Pj  £  J,  pl<p2’  where  4>(p>N,  [p-^Pj])  is  the 

fraction  of  the  iterates  {p , . . . ,G^~^(p) }  of  p  £  tpi,p2^’ 

Moreover,  the  set  =  {q  :  g(q)  >  0}  is  an  interval,  and 
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is  the  positive  limit  set  of  almost  all  initial 
point  pe  J. 

The  mapping  in  Figure  5.1 (a')  satisfies  the  hypotheses 

of  (S.5.3).  From  (S.5.3),  the  behavior  of  { Gn ( p ) }  is 

described  schematically  as  follows.  Choose  an  arbitrary 

pej.  As  n-*-°°,  Gn(p)  wanders  around  in  J  and  finally, 

(Gn(p)}  fills  an  interval  J^.  The  density  gCp)  of  points  in 

the  limit  n  -*■ 00  is  unique  and  independent  of  almost  all 

initial  points  pe  J,  and  is  mapped  to  the  same  g(q)  by  G. 

We  can  predict  the  behavior  of  the  system  in  the  limit  n  ■+• 00 

only  in  some  probabilistic  sense  although  the  system  is 

deterministic.  For  example,  we  know  that  the  probability 

for  the  point  to  be  in  [p^Pj]  in  the  limit  n  <=°  is 

.  Very  little  is  known  for  other  mappings  which  do 

not  satisfy  the  hypotheses  of  (S.5.3). 

The  mapping  in  Figure  5.1(c')  is  peculiar,  since  it 

has  a  countably  infinite  number  of  discontinuities.  If  we 

choose  T.  =  {X  :  x,  >  x^  ,  x.  =  x„  }  and  T  =  {X  :  x.  <  x7  , 

+  1  le  2  2e  -  1  le 

x2  =  x2e^’  t^ien  we  have  a  one-dimensional  mapping  in  Figure 
5.1(f)  and  its  corresponding  mapping  (ff)  for  the  simple 
model.  Here,  the  domain  of  mapping  is  an  unbounded  interval. 
We  assume  that  a  density  g(p)  of  p  exists  for  this  system, 
then  g(p)-»-0  as  p  -*■  0 .  We  know  that  the  magnitude  of  any 
physical  parameter  must  be  finite.  Here,  p  can  be  infinite, 
but  with  zero  probability.  Hence,  there  seems  to  be  no 
contradiction  to  the  requirement  as  a  physical  model.  Such 


J?2g(q)dq 

P1 
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an  interpretation  may  provide  new  types  of  physical  models. 
For,  the  magnitude  of  variables  have  been  conventionally 
classified  only  two  cases,  that  is,  bounded  or  diverging  to 
infinity,  and  no  concept  of  probability  has  been  involved 
for  such  a  simple  deterministic  model. 

5.3  Statistical  Properties  of  (3)  and  (6) 


Quasi  one-dimensional  mappings  obtained  experimen¬ 
tally  from  (3)  and  (6)  are  useful  to  understand  the  behavior 
of  (3)  and  (6).  But,  as  mentioned  in  the  previous  sections, 
it  is  difficult  to  study  their  statistical  properties 
analytically.  The  only  one-dimensional  mapping  whose  statis¬ 
tical  property  is  well-known  [2  2]  is  given  by 


G(p)  = 


2p ,  0  <  p  <  1/2, 

2 ( 1-p) ,  1/2  <  p  <  1. 


(5.1) 


It  is  known  that  if  pe  [0,1]  is  described  by  a  binary 
sequence  {b^},  the  map  G  acts  on  {b^}  like  a  shift  map  which 
shifts  terms  of  the  sequence  to  the  left  (Bernoulli  shifts). 
A  Bernoulli  shift  is  a  stochastic  system  (i.e.,  the  system 
behavior  is  completely  unpredictable  in  a  certain  sense 
[33]).  For  other  maps,  very  little  is  known.  Here,  we  try 
to  obtain  the  statistical  properties  of  (3)  and  (6)  numeri¬ 
cally. 

5.3.1  Mixing  Property 
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If  there  is  an  attractor  A  which  is  ergodic  and  y(X) 
is  a  defined  measure  on  A,  then  for  any  smooth  function  u, 

G4  /au(X» )dy(X' )  =  lim  (l/T)/Jdtu(X(t ,XQ) ) 

-►OO 

for  almost  all  initial  conditions  XQ .  The  system  has  a 
mixing  property,  if  for  any  smooth  functions  u  and  v, 

R(t)  A  Qv  -  /au(X(t,X’ ))v(X' )dy(X' ) 

=  uv  -  lim  (l/T)/Jdtu(X(T+t,Xn) )v(x(t,x  ) ) 

T-h»  u  u  u 

for  almost  all  X^  and  R(t)>0  as  t  -►<».  As  mentioned  in 
Chapter  1,  this  property  implies  that  the  time  evolution  of 
the  trajectories  is  highly  sensitive  to  the  initial 
condition,  which  is  an  important  property  of  the  turbulent 
solutions.  The  system  having  a  mixing  property  is  a 
stochastic  system  in  a  weaker  sense  than  a  Bernoulli  shift. 
We  examined  numerically  R(x)  with  u(X)  =  v(X)  ix^  The 
results  for  Figures  5.1(a)-(e)  are  shown  in  Figures  5.2(a)- 
(e),  respectively.  We  can  not  try  all  smooth  functions  u 
and  v,  but  the  results  in  Figure  5.2  suggest  that  the 
systems  have  the  mixing  properties  for  certain  values  of 
parameters . 

5.3.2  Unbounded  Solutions  of  System  (6) 

As  mentioned  earlier,  system  (6)  with  y^=0.2  has  an 
unbounded  solution.  If  its  first  return  mapping  as  a 
discrete-time  system  has  a  density  g(p)  of  p,  g(p)  -*■  0  as 
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p -*■  0 .  As  a  continuous-time  system,  system  (6)  must  satisfy 
the  physical  requirement  that  the  time-averaged  energy  is 
finite.  We  consider  the  ratio: 

PT<X)  =  measure{t  :  x^t)  >  x>  0  <  t  <  T}/T. 

In  Figure  5.3,  we  have  plotted  those  values  of  x  for  which 
PT  is  insensitive  to  the  initial  conditions.  As  T  increases, 
the  interval  of  such  values  of  x  also  increases  in  size. 
Hence,  it  is  likely  that  as  T  +  00,  P,p  approaches  P^  which  is 
independent  of  the  initial  conditions  and  goes  to  zero  as 

Ixl  - 

3 

If  exists  and  decays  more  rapidly  than  1/x  ,  and  if 

E]l  =  lim  (l/T)/JXl(t)2dt  =  /Jx2P00(x)dx 

is  satisfied,  then  from  the  physical  point  of  view,  we  can 
state  as  follows:  the  energy  of  wave  1  of  system  (6)  can  be 
arbitrarily  large  over  a  finite  time  interval,  but  its  time 
duration  is  so  short  that  the  time-averaged  energy  E-^  is 
finite . 

Equation  (6)  is  obtained  by  neglecting  higher  order 
terms  with  small  coupling  coefficients.  Hence,  if  fl X ( t ) fl  is 
sufficiently  large,  model  (6)  is  not  valid.  From  the  above 
argument,  the  time  duration  for  which  x-^Ct)  >  x  becomes 
smaller  in  the  ratio  P,p(x)  as  x  becomes  large.  Hence,  if 
the  range  for  the  approximation  is  large,  we  can  state  that 
model  (6)  is  valid  most  of  the  time. 


CHAPTER  6 


CONCLUSIONS 


6.1  Summary  of  Results 

In  this  work,  we  have  studied  the  chaotic  behaviors  of 
wave-wave  interacting  systems  involving  at  most  four  waves 
(Table  2.2).  These  systems  can  be  described  by  a  real  three- 
dimensional  equation  on  a  certain  set  in  the  state  space, 
where  the  phases  of  the  waves  are  locked  for  . 2  *  >  —  C 7 ’ )  (Sec. 
2.2).  We  have  found  that  there  remained  only  six  systems  of 
positive-energy  waves  (Table  3.2),  if  we  discard  the  cases 
which  have  the  following  property  :  The  corresponding  reduced 
equations  have  only  one  equilibrium  point  or  a  stable 
equilibrium  point,  or  they  have  an  open  set  RQ  in  the  state 
space  such  that  if  Xq  e  RQ,  X(t)  =  X(t,XQ)  e  RQ  for  all  t  >  0 
and  diverges  to  infinity  as  t  +  ».  It  is  very  likely  that 
the  reduced  equations  (D  —  (4) ,  (6)  have  chaotic  solutions  for 
almost  all  initial  conditions  in  the  state  space.  System 
(3')  and  (6')  have  only  linear  damping  terms,  and  the 
instability  of  their  equilibrium  points  is  due  to  the 
external  wave.  The  numerical  experiment  was  done  only  for 
these  two  cases. 

We  obtained  the  conditions  (Table  2.7  and  (S.2.3))  for 
asymptotic  phase  locking  of  (2')-(7')  as  t-*-°°,  which  depend 
on  the  sizes  of  the  attractors  for  the  trajectories.  But, 
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the  conditions  are  too  strong  for  determining  the  occurrence 
of  asymptotic  phase  locking  of  (3')  and  (6')  as  t  -*■ 00  for  the 
values  of  parameters  at  which  chaotic  solutions  may  appear. 
The  numerical  experiments  show  that  system  (6')  has  a  strong 
tendency  of  asymptotic  phase  locking  as  t+®.  In  (3*),  for 
certain  values  of  parameters,  the  phases  are  nearly  locked 
most  of  the  time,  but  not  asymptotically  locked  as  t -*■“>. 

Periodic  and  pseudo-chaotic  solutions  of  system  (3) 
and  especially  (6)  were  studied  numerically  in  detail.  For 
certain  values  of  parameters,  the  first  return  mappings  of 
the  trajectories  have  properties  which  are  similar  to  those 
of  the  well-known  one-dimensional  mappings  having  parabolic 
graphs.  The  transitions  of  attractors  are  schematically 
explained  by  using  the  simple  models  in  Figures  4.7  and  4.8. 

Finally,  we  considered  the  statistical  properties  of 
the  first  return  mappings  of  the  trajectories  of  (3)  and  (6) 
using  the  known  theorems  for  simple  one-dimensional  discrete 
models.  Numerical  results  suggest  that  system  (3)  and  (6) 
have  a  mixing  property  for  certain  values  of  parameters, 
which  is  an  important  property  of  the  turbulent  states. 

6.2  Remarks  on  Further  Research 

We  have  shown  that  some  wave-wave  interacting  systems 
with  a  few  modes  have  chaotic  behaviors.  As  mentioned 
earlier,  no  additional  assumption  for  randomness  was 
necessary  in  these  systems.  In  this  new  approach  to  plasma 
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turbulence,  there  are  still  many  unresolved  basic  questions. 

One  of  the  questions  is  how  to  relate  the  statistical 
properties  of  the  simple  systems  with  a  few  modes  to  the 
turbulent  state  in  plasmas.  In  actual  plasmas,  there  are  a 
large  number  of  coupled  modes.  If  a  plasma  in  a  turbulent 
state  consists  of  systems  which  are  weakly  coupled  with  each 
other  and  each  of  which  is  a  few-mode  interaction  having  a 
chaotic  motion,  we  can  regard  each  system  as  a  quasi¬ 
particle,  since  the  phase  of  the  system  varies  randomly. 
Then,  the  whole  system  may  be  described  by  a  kinetic 
equation  of  quasi-particles  [11].  If  many  waves  are 
strongly  coupled,  we  cannot  use  directly  the  results  for  a 
few-mode  system  to  analyze  the  whole  system.  It  has  been 
shown  that  a  multi-mode  system  has  behavior  similar  to  that 
of  a  few-mode  system  [10].  But,  very  little  is  known  about 
general  multi-mode  systems. 

There  are  also  unresolved  questions  about  the  simple 
systems  studied  in  this  work.  First  of  all,  we  still  do  not 
know  exactly  under  what  conditions  the  reduced  equations 
describe  the  asymptotic  behaviors  of  the  original  systems. 
Furthermore,  the  relationship  between  the  reduced  and 
original  systems  is  not  known  if  the  phases  are  not  locked 
asymptotically  as  t  -*<*>.  None  of  the  statistical  properties 
mentioned  in  Chapter  4  are  proved  mathematically.  The  power 
spectrum  of  the  variables  is  one  of  the  important  quantities 
from  the  physical  standpoint.  But,  at  the  present  time,  it 
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is  too  difficult  to  obtain  the  power  spectrum  of  the  chaotic 
solutions  of  such  systems  analytically.  Thus,  further 
studies  must  be  done  before  the  turbulent  behavior  of 
plasmas  can  be  analyzed  via  the  new  approach  considered  in 
this  work. 
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APPENDIX 


PIKOVSKII-RABINOVICH-TRARHTENGERTS ’  MODEL 

Let  us  assume  that  all  the  y^'s,  v^'s  and  v|'s  are 
positive.  By  using  the  substitution: 

A^  +  iA^,  A2+A2J  A^  +  A^,  vAg-*-h,  (A.l) 

we  obtain  the  following  set  of  equations 

A1  =  -  Y-jA^  -  A2A3  +  hA*  , 

A2  =  -  Y2A2  +  AiA*  +  ’  (A. 2) 

A3  ~  —  Y 3A3  *  AqA2 

which  is  identical  to  the  Pikovskii-Rabinovich-Trarhtengerts ' 
model  [12],  Let 

5  =  Im(hA*A*),  n  =  ImCAjA^A*).  (A. 3) 

Then, 

(d/dt)(5  +  n)  =  -  (Yj  +  Y2  H  -  +  Y2  +  Y3^n’  (A. 4) 

According  to  Ref. [12],  the  trajectories  S =  (£(t) ,n(t) )  enter 
the  sector  Z  between  the  straight  lines  L^  and  L2 : 

L1  =  {(£,n)  :  S+n  =  0), 

(A. 5) 

L2  =  {(£,ri)  :  (Y1+Y2)5  +  (Y1+Y2+Y3>n  =  0}.  , 

It  is  true  that  the  trajectories  E=  (£(t),n(t))  on  L^  enter 
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Z,  since 


(d/dt)(5*n)|(?>n)e  Li=‘-T3n, 


CA.6) 


which  is  positive  (resp.,  negative)  when  n  is  negative 
(resp.,  positive).  However,  the  trajectories  S = (£(t) ,n(t) ) 
on  L>2  do  not  always  enter  Z.  For, 

(d/dt){(Y1+T2H+(Y1+Y2'HY3)r|^  =  (V1+Y2+Y3KY1+Y2H?+n) 

-Y3(2Y1+2Y2+Y3){<Y1+Y2K+<Y1+Y2+Y3)h} 

+hY3{Im(A12A*)-Im(A22A3)}  =  Q,  (A.7) 


and  Q  is  easily  shown  to  be  positive  (resp.,  negative)  on 
for  £>0  (resp.,  £<0).  Without  loss  of  generality,  we  set 
h>0  since  the  substitution  h-*-h,  A^  -*■  -A^  and  A3  -*>  -A3  leads 
to  the  same  equation  (A. 2).  On  Lj,  we  can  choose  any  values 
of  |A1 1  and  J A^  J  independently  of  r3  and  Qj's,  since 

r1r2{-h(Y1+Y2)sin(01+92)+r3(Y1+Y2+Y3)sin(ei-e2-03)}  =  0  (A.  8) 


on  L2 .  We  can  also  choose  the  S^'s  such  that  sin(20^-02)  i  0 
and  sin(282+03)  i  0,  since  (A. 8)  has  two  degrees  of  freedom 
for  the  0j's.  Now,  we  assume  that  £  >  0,  then  n  <  0  and 
S+n  >  0  on  L2 .  On  L2 ,  we  can  choose  only  0j's  such  that 


sin(20^-©3)  <  0  and  sin(202+03)  <  0,  or 
sin( 20^-©3 )  > 0  and  sin(  2 0 2  +  0 ^ )  >  0  . 


(A. 9) 


Hence , 
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(A. 10) 


Q =  (Y1+Y2)(Y1+Y2+Y3)rlr2lhsin(9l+e2)+r3sin(ei“92“e3) I 

+hY3{+r12r3|sin(201-e3) |+r22r3jsin(262+03) | }. 

Thus,  for  sufficiently  large  |r2|/|r^|,  Q  is  positive  on  L2 . 
Similarly,  Q  is  negative  on  L>2  for  £  <  0  for  sufficiently 
large  |r^|/|r2(.  Hence,  trajectory  (£(t),n(t))  may  leave 
sector  Z. 


